T 1 044F 


Applied Dynamics Research Corporation 
Final Report 

Investigation of Empirical Damping Laws 
for the Space Shuttle ■ 


Contract NAS8- 2851 3 14 March 1973 


Prepared by 

Edward L. Bernstein 
Research Specialist in 
Structural Dynamics 




Approved by 






C. S. Chang 

President 

ADRC 




FOREWORD 


This report presents the work performed by Applied Dynamics 
Research Corporation under Contract NAS8-28513 for the Marshall Space 
Flight Center of NASA. Mr. Larry Kiefling, S&E-AERO-DDS is the 


technical monitor. 



SUMMARY 


An analysis of dynamic test data from vibration testing of a num- 
ber of aerospace vehicles is made to develop an empirical structural 
damping law. A systematic attempt is made to fit dissipated energy/ 
cycle to combinations of all dynamic variables. The best-fit laws for 
bending, torsion and longitudinal motion are given, with error bounds. 

A discussion and estimate are made of error sources. Programs are 
developed for predicting equivalent linear structural damping coefficient 
and finding the response of nonlinearly damped structures. 
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Section 1 


INTRODUCTION 

A number of full scale dynamic tests on aerospace vehicles have 
been made which produced enough data for the calculation of damping 
properties. Some preliminary analyses of the data (Ref. 1 and 2) indi- 
cated a correlation between dissipated energy per cycle and the total 
vibratory energy of the vehicle. Based on the success of these analyses, 
and a number of suggestions proposed (e. g. , Ref. 3) on other parameters 
which might influence damping, a comprehensive study was made to 
determine an empirical damping law. 

Dissipated energy/cycle was selected as the physical quantity 
which provides the best measure of damping. Damping in a complicated 
structure is a combination of a number of types of phenomena, and is 
often nonlinear. When an elastic structure is excited periodically, a 
steady-state situation is reached when the rate of energy dissipation 
through the damping mechanisms becomes equal to the rate of work put 
into the system by the excitation force. The rate of work can be calcu- 
lated independently of the type of damping, and is therefore chosen as the 
logical dependent variable for an investigation of structural damping. 

The work done, per cycle of vibration, can be calculated as 
follows: Let 
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f(t) = F sin co t 


denote the exciting force, and 

X E = X E s ^- n ® ) 

be the displacement at the point of application of the force. Here, co is 
the frequency of the excitation force, and 0 is the phase angle of 

the response, relative to the phase of the excitation. 

The work per cycle, W, is then by definition 


W = 


2y 

a 

•fn 


f(t) x E dt 


= 7rF X E sin 0 

At resonance, 0 = tt/ 2. Thus, 

W = ttFX_ 

£ 

which, from the previous discussion may be set equal to the energy 


dissipated per cycle, D: 

D = ttFX e ( 2 ) 

An alternate method of calculating D can also be used, based on 
the value of the equivalent linear damping coefficient f • The damping 
at a particular modal amplitude is described by this quantity, although 
it will not be valid at other amplitudes if the structure is nonlinear. The 
dissipated energy is then given by 

D = 4 7rf T ( 3 ) 

where T is the maximum kinetic energy. ; 


2 



Dissipated energy was calculated from both formulas when suffi- 
cient information was available. A comparison of the two values provided 
a measure of the accuracy of some of the basic data used. 

The data available for this study included the results of a number 
of full scale dynamic tests on space vehicles, as well as the results of 
tests on a Boeing 747 prototype. These tests were performed to deter- 
mine natural modes and frequencies of these vehicles. Most tests of 
this sort were not run with the intent of obtaining nonlinear damping data. 
Therefore, as far as this program is concerned, useable damping data 
were relatively limited. In addition, the accuracy of some of the damping 
data produced is poorer than the accuracy of the frequency determination 
alone, which generally was the main purpose of the experimentation. 

Most of the tests were run on Saturn vehicles, and an evolution through 
the years of increased accuracy of methods for determining the dynamic 
parameters is seen. 

A considerable attempt was made to locate data on dynamic 
analyses of solid-fueled rockets which might be most applicable to the 
space shuttle. This attempt, however, was unsuccessful. As far as 
could be determined, the data were unavailable. 

The data available for this study comprise the results of the 
tests given in Table 1. Altogether, there are 15 tests representing seven 
structures. Five of the structures are configurations or stages of the 
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Saturn V. The Saturn I tests were run in 1964. Damping constants were 
available for some of the modes from ringout decay tests. This data was 
not considered too reliable due to the beating phenomenon which frequently 
occurred. Generalized mass was obtained by measuring the bandwidth of 
the half-power point on resonance peaks of the amplitude /frequency curves. 

The S-IVB tests were run in 1965-66. Generalized mass was found 
from the mode shape and mass distribution and also from a complex curve- 
fit for the amplitude of the response curve. Damping values were obtained 
from bandwidth and phase angle methods, and, for some tests, from a 
complex curve -fit technique. 

The Saturn V tests were performed in 1966-67, and represent a 
well-documented and more accurate effort. Generalized masses and 
damping constants were found from complex curve -fit techniques. 

The Boeing 747 tests were performed in 1969. Dynamic proper- 
ties were determined from complex admittance (acceleration divided by 
force) plots in polar coordinates, following the method of Kennedy and 
Pancu (Ref. 4). 

The first indication of the existence of an empirical damping law 
arose from analysis of the Saturn I data (Ref. 1). In this analysis, 118 
data points were used. The reason for elimination of the remaining data 
is not known. The missing data mostly consist of entire tests in the 
series. Very likely, at the time of the analysis, facts were available 
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indicating the unreliability of these particular tests. The conclusion 
reported in Ref. 1 was that a law of the form 

D = 0. 313 T 0,8 (4) 

[dimensions in SI (International System) units] predicted structural damp- 
ing. It was reported that 83% of the data lay within a _+ 2 db band of this 
law. The data used was reexamined, using the methods of the present 
study. It was found that the law was not quite so good as reported. The 
best law for the S-I data was found to be 
D = 0. 350 T 0 * 72 

with 75% of the data within a + 2 db band. This result still is encouraging 
for the existence of a structural damping law. 

Based on the S-I analysis, Kiefling and Pack (Ref. 2) analyzed 
samples of data points from the Saturn V and the S-IVB tests. From 
these tests, data on torsional and longitudinal vibrations were also avail- 
able. The results reported in Ref. 2 indicated that the law (4) predicted 
damping in bending quite well for these structures. As a matter of fact, 
the plotted results in Ref. 2 indicate that 89% of the data fall within a 
+ 2 db band. 

For torsional vibrations a law of the form 
D = 0. 34 T°* 80 

is reported in Ref. 2, with 67% of the data within a + 2 db band. The 
results for longitudinal vibrations are reported as inconclusive. 
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In Ref. 5, Riley took vibration data from a Boeing 747 airplane 
and applied the Chang law (Eq. 4). He found that this law correctly pre- 
dicted damping to within + 2 db for 75% of the data. 

While the accuracy of the law, based on these preliminary studies, 
is not all that might be hoped for (+ 2 db represents a + 50% and - 37% 
error in D), yet it represents a remarkable correlation. A wide variety 
of vehicles are included (single tank, multi-tank and winged) and a five- 
decade range of energy levels. Taking into consideration the fact that 
dissipated energy is an overall quantity, depending on an unknown variety 
and number of damping mechanisms acting within each structure, it 
therefore appears that an investigation to uncover a damping law is 
worthwhile. 

The studies reported above indicate a dependence of dissipated 
energy on total vibration energy. Other suggestions have been made that 
frequency to may have some influence. In Ref. 3, it is reported that 
the ratio of amplitude x to some characteristic vehicle length L gives 
good correlation for certain cases. 
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Section 2 


SCOPE AND METHOD OF THE INVESTIGATION 

In order not to exclude any possible empirical laws, dependency 

on almost all variables available from dynamic testing was investigated. 

In addition, all possible combinations of these variables were tried. The 

variables used, together with their units in the SI system, are 

T -- kinetic energy (Newton- meters) 

X -- mode shape amplitude (meters) 
co — frequency (rad/sec) 
m -- generalized mass (kilograms) 

X/L -- amplitude /length 

These quantities are not all independent, since the kinetic energy 
is given by 

T = (l/2)mw Z X 2 (5) 

Correlations were sought for three cases: transverse bending, 
torsion, and longitudinal motion. The length L, in the last of the inde- 
pendent variables, was chosen to be the length of the structure in the 
case of bending and longitudinal motion, and maximum radius for the 
case of torsion. 

Equation (4) for the empirical damping law is not independent of 
the dimensional units of the’ variables. - Some effort was made to find a 
satisfactory nondimensional formulation. The effort, however, was 
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unsuccessful. Therefore, for the sake of uniformity, and to conform 
with current standards, all units were converted to the SI system. 

The energy levels in these tests covered a five-decade range. 
Based on the preliminary analyses, a law was sought of the form 

D = CP 1 ai P 2 a2 P 3 a3 ( 6 ) 

where C, a^, a 2 , a^ are constants to be determined, and Pj, P 2 , P3 
are independent parameters. 

Taking the natural logarithm of both sides, we have 

d = In D = c + ajPj + a 2 p 2 + a 3 p 3 ( 7 ) 

where 

c = In C 

Pp p£> P3 “ ln ( p i» p 2’ p 3) 

The least- squares criterion was used to find the constants. Let 
the data points determined from the tests be denoted by a superscripted 
bar. Then the error, E, is given by the expression 

n 

E = ^ (c + a^j. + a 2 p 2 . + a 3 p 3 . - d.) 2 (8) 

• = 1 

where 

d-i = ln Dj 

Pli* P2P P3i = ln < P li' P 2i’ P 3i> 
and the summation is over N data points. 
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We note that on a log-log scale the sought-after empirical law is 
a straight line. The optimal values of C, aj, a 2 , a 3 are found from the 
equations for minimizing E: 

3E _ 9E _ 9 E _ 9 E _ q ^ g j 

9c 9 el ^ 9 ^2 9 a 3 

Upon substitution of (8) into (9), we have 

E (c + a lPu + a 2 ? 2i + a 3^3i - 3 i> = 0 

E Pli <C + a lPli + a 2 ? 2i + a 3P3i - V = 0 

(10) 

^P 2i (c + a 1 p 1 i+a 2 P2i + a 3 P 3i -^)= 0 

^p 3i (c + a^. + a 2 p 2 . + a 3 p 3 . -3^=0 

The simultaneous solution of these equations yields the best-fit values 
of C , a^ , » a^. 

The following formulas were used in searching for a best fit: 

D = C T a l 
D = C X n ai 
D = C co a l 
D = Cm a l 
D = Cm a lX a 2 

D = Cm a lo; a 2 (11) 
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D = CX a lco a 2 
D = CT a l(X/L) a 2 
D = Cm a lco a 2x a 3 
D = Cm a lco a 2(X/L) a 3 
D = Cj + C 2 T + C 3 T 2 

A program was written for the digital computer (an XDS 930) which 
performed the least- squares analysis and calculated the error of the final 
result. The program performed several functions. Data in a number of 
dimensional units were read into the program and converted into SI units. 
Dissipation energy and maximum kinetic. energy were calculated. Each of 
the formulas (11) was used with the least- squares criteria (10), and best- 
fit values for the unknown constants found. Finally, the error was calcu- 
lated by finding the percentage of data points which fell within + 1 db, 

+ 2 db, + 3 db, and + 4 db bands of the empirical law. A listing of the 
program is given in Appendix B. 

As mentioned previously, there are two ways to calculate D. 
Equation (2) is the preferred way, since it is a direct calculation involving 
measured dynamic parameters. However, measurements of the ampli- 
tude at the location of the excitation force may be inaccurate due to the 
small magnitude there. For .those tests where the damping constant was 
available, the dissipated energy was calculated according to both equa- 
tions (2) and (3). For each test, the computer program was used to find 
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the best-fit law of the form 
D = CT n 

and the scatter calculated (percentage of data points varying by + 1 db, 

+ 2 db, + 3 db, + 4 db). In this way, an empirical measure of the accuracy 
of both formulations could be gauged. This subject will be more fully dis- 
cussed in a subsequent section on error analysis. The conclusion reached 
was that the calculation of D from equation (3) was generally more accu- 
rate. This was markedly noticeable for the S-IVB tests, and to a much 
lesser extent the case for the S-V tests. 

In the Saturn V series of tests, dynamic parameters were meas- 
ured in two types of tests. The primary test, conducted to find the dynamic 
parameters, was a frequency sweep test. An additional test run was a 
force level test, mostly for the purpose of checking linearity of response 
over a range of excitation forces. The results of this test, however, did 
supply sufficient information to calculate the dynamic parameters. In 
the force level tests, responses were found at a large number of vehicle 
locations for three force levels. The maximum level was approximately 
the same as in the frequency sweep test. For this magnitude of force, 
the results of the two types of test were compared for consistency. 

The results of the force level test were potentially quite valuable, 
since it presented an opportunity to determine damping in a situation 
where all but one parameter remained constant. However, for the 


11 



torsional and longitudinal tests, there was poor correlation between the 
force level tests and the frequency sweep tests for the same excitation 
force. In addition, the scatter was quite large. Therefore, the force 
level tests for torsion and longitudinal motion were not used in this 
analysis. 
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Section 3 


RESULTS 

The results are given in Tables 2-5, which cover bending, torsion 
and longitudinal motion. Each table gives the damping laws which fit the 
data best, for each of the formulations tried from equations (11). The 
scatter is given by the percentage of data points falling within _+ 1 db, 

+ 2 db, + 3 db, and 4 db bands around the basic law. Figures 1-4 show 
the data as a function of kinetic energy, with the line representing the 
damping law. 

The laws for bending represent an analysis of 330 data points, 
much more than the results for torsion or longitudinal motion. For this 
reason, it is felt the results are more accurate and likely to be of general 
application. The conclusion derived from Table 2 is that a law of the 
form 

D = 0. 286 T°‘ 746 (12) 

describes the damping energy of a large class of aerospace structures. 
Figure 1 illustrates the fit of this law to the data. 

Of the 330 bending data points, 196 or 60% are associated with 
Test No. 1. Therefore, the damping law (12) may be biased by the 
inclusion of this data. For this reason, an alternate damping law for 
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bending was derived, based on tests 2, 5, 8, 11, 13 and 15. The law 
derived was 

D = 0. 153 T°' 893 (13) 

with 71% of the data falling within a + 2 db band. The results of the in- 
vestigation are given in Table 3. Figure 2 illustrates this law. 

It is interesting to note that when a damping law was sought of the 

form 

D = Cm a lw a 2x a 3 

the law which fits the data best produced the approximate relationship 
a ^ ft? 2a^ 

•y 

That is, the quantity mX^ influences the damping. A glance at Tables 2, 3, 
4 and 5 reveals this is true for bending, torsion and longitudinal motion. 
Although a law of this form gave the best fit to the data, it was not 
significantly better than the simpler law (12)- -not enough to recommend 
its use. 

For torsional vibration, the best law was found to be 
D = 0.095m 1 - 01 co 2 - 69 X 2 - 33 

with 81% of the data falling within a _+ 2 db band. The law 

D = 0. 101T 1, 279 (14) 

had a 63% accuracy for + 2 db. While the former law shows a better fit, 
it should be borne in mind that the total number of data points is only 26, 
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so that a small uncertainty in the five points which lie just outside the 
+_ 2 db band (see Fig. 3) would change 63% accuracy to 81% accuracy. 

The law for torsion is less reliable than the bending law, since 
it represents an analysis of only 26 points, all relating to configurations 
of the Saturn V. It is noteworthy that the law (14) for torsion has an 
exponent greater than unity, implying damping ratio increasing with 
amplitude. This is the reverse of the situation for bending. 

For longitudinal vibration, the best law was found to be 

D = 0. 0002m 1 ' 259 co 3, 315 X 2-405 (15) 

with 77% of the data falling within a + 2 db band. The law 

D = 0. 057T 1 ’ 1 04 (16) 

had a 67% accuracy for + 2 db. Here, the number of data points is 39, 
all representing the Saturn V with and without the SI-C stage. Figure 4 
reveals that a small uncertainty in the. data could give as much as 77% 
accuracy for the law (16) 

Again, the law for longitudinal motion has an exponent greater 
than one. 
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Section 4 


ERROR ANALYSIS 

One cannot help noticing that the empirical laws shown in Figures 
2-5 are subject to a considerable degree of scatter. The natural question 
to ask is how much accuracy can be inherently expected from the data 
generated. The tests were not primarily designed for the purpose of 
measuring dissipated energy, and a careful look at the origins and accu- 
racy of the dynamic parameters is therefore in order. 

The quantities plotted in Figures 2-5 are derived from the formulas 
T = (l/2)mo; 2 X 2 (5) 

and 

D = f T (3) 

or 

D = ttFX e (2) 

First, it is necessary to derive expressions for the error in the 
derived quantities T and D propagated from the independently determined 
quantities m, co , X, X E> f , F. Of course, some of these quantities 
themselves are derived from other independently determined variables. 

In general, if there is a continuously differentiable function 
U = U(x, y, z) 
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and x, y, z are replaced by their approximate values 

x = x + e x 

y = y + e Y 

z = z + e z 

where e x , Cy, e z are the errors in x, y, z, then the error in U is e 
given by (Ref. 17) 

au , au , au ^ 

e u = ax Cx + ay e y + az z 


u 


(17) 


Applying equation (17) to formulas (2), (3) and (5), we have the 
following expressions for relative error: 


e T e m 


+ 2 — — + 2 — - 
m w x 


(18) 


and 


D 

D 


+ 


T 

T 


(19) 


or 


_ Jx. + e x E 


D F X E 

Thus, if some measure can be assigned to the errors of the 
independent variables (independent with respect to T and D), then the 
errors in T and D are known. , 

Of the several series of dynamic tests, the Boeing tests on the 
Saturn V (tests nos. 2-10) are the best documented and, at least on the 
basis of currently available information, the most accurate. The 
accuracy of these tests is discussed below. 


( 20 ) 
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For many of the tests, the energy dissipated/ cycle could be 
calculated from either formula (2) or (3) independently. A comparison 
of the values calculated is shown in Table 6, as a measure of accuracy. 

It is seen that agreement is poor for many cases. In the bending tests, 
however, cases 2, 5 and 8 (the Boeing tests) show relatively good accu- 
racy. The discussion of probable error in the recorded values of dissi- 
pated energy and kinetic energy will begin with the Saturn V tests in 
bending. 

For the Saturn V bending tests, References (18) and (19) provide 
a source for estimation of errors. Errors accumulate from sensor cali- 
bration, instrument resolution, hysteretic effects in vehicle response, 
noise and scalar errors in the data acquisition system, and round-off and 
approximation errors in the data reduction methods. 

The Boeing tests incorporated a number of accuracy-correction 
techniques. Calibration curves for the amplitude and phase of each sensor 
were determined and used in the data reduction system to modify the final 
output. In addition, on-site calibration was periodically performed to 
check drift and malfunction. A Fourier analysis was performed on the 
response to remove noise and harmonic content. The output data (in the 
form of transfer functions for each sensor) is curve -fitted to yield the 
dynamic parameters. 
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The measurement of force was subject to two sources of error. 

It is not possible to assign a value to the error, which could have been 
important in the evaluation of dissipated energy for those cases where 
formula (2) was used. 

Misalignment of the load cell can cause significant error in the 
force reading. Misalignment will occur to some extent because of the 
bending curvature of the structure as it undergoes deformation. A second 
cause of misalignment is due to the out-of -plane motion of the structure. 
Since the excitation was not applied in a principal plane, the vehicle re- 
sponse was at an angle with respect to the force. Therefore, the axis of 
the load cell may suffer mis orientation and show only a component of the 
applied force. This component will generally be quite small. However, 
the out-of-plane component can cause a moment acting on the load cell, 
depending on its displacement from the shaker, which is a potentially 
more serious cause of error. 

Another source of error is due to the mass of the shaker armature 
causing an effective inertial force showing up in the load cell readings. 
The total test system should properly be considered to be the test vehicle 
plus the shaker armature. A mass compensation procedure should be 
performed to eliminate the effect of the additional mass. No indication 
was found in the documentation that this procedure was used. An example 
of the possible error from this cause is as follows. 
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For the Boeing tests, a 20,000 lb. (88,965 N) thruster system was 
used. A system of this magnitude uses about a 350 lb. (159 kg. ) arma- 
ture. For a frequency of 8. 69 Hz, and a thruster travel of 5. 35x10" m. , 
the corresponding inertial force is 253 N. This compares to a measured 
force of 5791 N, or a 4.4% error. It is therefore possible, especially at 
the higher frequencies, that lack of mass compensation could be a signifi- 
cant source of error in the force measurements. 

One of the principal objects of the dynamic tests was to determine 
modal natural frequencies. It is to be expected, therefore, that the accu- 
racy associated with this quantity is relatively good. There are several 
sources of error: 1. Numerical error associated with the curve-fit 

procedure. In Reference (18) an error analysis of the curve-fit technique 
is performed. For the one example given, the error in amplitude varies 
from 0. 1-0. 7% over the frequency range. 2. The possibility of a double 
peak in the vicinity of the resonant frequency. This effect is due to the 
fact that the excitation is not applied in a principal plane, and therefore 
the response is a summation of motion in two planes, each of which may 
have a slightly different peak due to asymmetry of the vehicle. In this 
case, as pointed out in Reference (19), the curve-fit method will not 
separate modes closer than 1% in frequency. 

From the above discussion, a value of 1. 5% is assigned, in a 
purely estimatory manner, as a typical error to associate with 
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experimental frequency values. Sensor inaccuracy may have little effect 
on frequency, since only the location of an amplitude peak is sought, and 
not its value. 

Errors associated with X (response amplitude) are influenced by 
sensor error and data acquisition and reduction error. The sensor error 
depends on the response magnitude. According to Reference (19), sensor 
error averaged 5% for measurements at 1% of full scale, while for full 
scale measurements, they averaged 1%. 

A check of the accuracy for determination of response amplitude 
X comes from comparison of the results of the frequency sweep tests 
with the force level tests. The force level tests were conducted as a 
check on the linearity of response as a function of force. The force levels 
of the frequency sweep tests were approximately reproduced, as well as 
two lower levels, and responses at these values recorded. Obviously, 
we expect agreement between the two tests at the same level. Compari- 
son of the results of the two tests reveals that the average agreement 
for X is within 5%. 

The determination of generalized mass is subject to the greatest 
inaccuracy of all the quantities making up the kinetic energy. This quan- 
tity is derived from the curve -fit technique. According to Reference (19), 
comparison of generalized mass values gave agreement to within 10%. 
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The modal damping values, likewise determined from curve fitting, 
were found in Reference (19) to give agreement within 5%. 

The above error figures maybe used to give an order -of -magnitude 
estimate of the error involved in the dynamic parameters of the Boeing 
tests. From equations (18) and (19) it follows that the error in determina- 
tion of kinetic energy is 

-lx- = 0. 1 + 2(0. 015) + 2(0. 05) = 0.23 = 23% 

and 

= o. 05 + 0. 23 = 0. 28 = 28% 

D 

This may be compared with the variation due to the + 2db bands in 
Figures 1 and 2 . Since the definition of a 2 db difference between two 
quantities Tj and T 2 is 

T, 

2 db-» log qr- = °-2 

then it follows that 2 db corresponds to an error of +59%, -37% in T. 

It is interesting to note that if we assume a Gaussian distribution 
of error, and that the nominal errors estimated represent a variation 
from the mean error of one standard deviation, then 68% of the data would 
lie within these error bounds. This figure correlates well with the 71% 


: 
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of data falling within the + 2 db bands of Figure 1 and 60% of data in 
Figure 2. 

The Boeing tests include 38% of the data points for all the dynamic 
bending tests, or 71% of all data points excluding the Saturn I data. 

The other data available on the remaining dynamic tests is gener- 
ally lacking in enough information to estimate the accuracy. This data 
falls into three categories. Test No. 15 was a ground vibration test of 
a Boeing 747 airplane. A frequency sweep was made, and frequency 
determined from polar plots of complex admittance (Ref. 5). From this 
plot, generalized mass and damping coefficients are also determined. 
From the rather sketchy details available concerning the conduct of the 
tests, no separate determination of accuracy can be made. 

Tests 11 and 13 were dynamic tests of the Saturn upper stages, 
performed by the Chrysler Corporation. The resonant frequencies were 
determined from inspection of the response and phase angle plots. Damp 
ing coefficients were found from an average of frequency bandwidth and 
phase angle methods. Generalized mass was found from integration of 
mass and mode shapes, using 70 equally spaced points along the vehicle. 

Not enough information was available to accurately assess the 
error bounds. Inspection of the response and phase angle plots indicates 
an up to 2% possible uncertainty in determining the resonance point. The 
generalized mass was found from the integration method mentioned and 
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compared with values found from a complex curve-fit for the amplitude 
of the response curve at the nose of the vehicle. Agreement between the 
two values was very poor. The integration method seemed to produce 
more consistent values and was the basis for the values used in this 
report. An estimated 15% error is associated with these values. The 
Chrysler tests account for 1 1% of all data points when the Saturn I is 
included, and 21% without the Saturn I data. 

Finally, the Saturn I tests, which formed the original basis for 
the Chang law, consist of 196 data points, or 59% of all data points con- 
sidered. Apparently, frequencies were found from inspection of response 
and phase angle plots. Generalized mass was found by integration of the 
mode shapes with a lumped-mass model. Damping constants were not 
used for determining dissipated energy, but rather equation (3), using 
force and excitation point response measurements. Insufficient informa- 


tion is available for error estimation. 



Section 5 


CALCULATION OF EQUIVALENT LINEAR 
DAMPING COEFFICIENT 

The results of the preceding section indicate that structural damp- 
ing in bending can be predicted by equation (12) or (13). With less cer- 
tainty, equations (14) and (16) may also be useful for predicting damping 
for torsion and longitudinal motion. Damping, at a particular amplitude, 
may be described by an equivalent linear damping coefficient f , given by 


From equation (17), the coefficient for bending is 
? = 0. 0228 T“°* 254 

For the range of T covered in the tests (see Fig. 1), this equation gives 
the variation of f to be 

0. 07 > f > 0. 004, for 0.01 N-m< T < lOOON-m 
For torsion, from equation (14), the coefficient is 
f = 0. 00804 T°* 279 

and gives a variation over the range of T (Fig. 2) of 

0. 004 < f < 0. 029, for 0. 1 N-m < T < 100 N-m 
Finally, for longitudinal motion, equation ( 1 6 ) gives 
f = 0. 00454 T°* 104 
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with f varying (Fig. 3) as 

0. 006 < f < 0. 009, for 10 N-m < T < 1000 N-m 
A small digital computer program was written to calculate for a 
range of amplitudes of kinetic energy. Therefore, the results of modal 
analyses of the space shuttle, or other large aerospace vehicles, can be 
used to produce predicted values of f for further application in response 
studies. A listing of the program is given in Appendix B. 
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Section 6 


RESPONSE OF NON LINEARLY DAMPED STRUCTURES 

The results for the equivalent linear damping coefficient, found in 
the preceding section, can be used to find the complete dynamic charac- 
teristics of a complex structure. We let the response at a point x of the 
structure by U(x,t). Then U may be expanded into a series of modal 
functions of the form 

U(x,t) = £q n (t)$ n (x) (22) 

where 

4> n (x) is a mode shape function, and 

q n is the generalized coordinate corresponding to 

The modal equation of motion for the n-th mode is 

% + W n q n + F * /m n = Q n (x) (23) 

where 

0 J n is the natural frequency, 
m n is the generalized mass, 

F n is a nonlinear function including damping, and 

Q n (t) is the generalized force, written, for periodic excitation, as 

Q n = Qn cos 00 t 

The fundamental response is 
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qnft) - q n ( w ) cos (w n t - 6-n ) 

where 6 n is the phase angle. 

F n is expanded into the first terms of a Fourier Series as 

F n =• I x cos (w n t - Oj + I 2 sin (w n t - $ n ) 

Dn 

= Ij cos (oo n t - e n ) - _ sin (co n t - 0 n ) 

nq n 

_ 2 %. 


(24) 


i 1 f __ 

where I = — / F„ 

n ir J n 


sin (o> n t - 0 n ) a> n dt 


and D n is the energy dissipated/cycle. Substitution of equation (24) into 
equation (23) gives the two equations 


q Ju„ 2 - oj 2 ) +IL = Q„ cos 6 


A n n 


m n 

m n 


n 


(25) 


D r 


— = Q sin 6 n 
^n m n 


The peak kinetic energy/cycle, T n , may be written in terms of q n as 


n 


- 1 2 n 2 

— — r.i a 


2 w q n m n 


Substitution of this equation into equations (25) and recombining gives 


q n ~ Q n [ (w n " ^ + + (D n w 2 /2 tt T n ) 2 J 


- 1/2 


2^ . 2 2 


tan<? n = w D n /2uT n (co n - w + ^ / m n q n ) 


28 



If the natural frequency co n does not vary significantly with response 
amplitude, I, = 0. Then q n and 0 n can be found by the simultaneous solu- 
tion of the nonlinear equations (25). • 

A digital computer program was written to accomplish this. The 
program requires as input the dynamic modal parameters w n , m n and the 
amplitude and frequency of the applied force. It then solves equations (25) 
iteratively through the use of Newton's method. The resulting solutions 
are then used in equation (22) to find responses for arbitrary forcing 
functions. A listing of the program is given in Appendix B. 
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Section 7 


CONCLUSIONS 

This study represents a comprehensive analysis of all available 
dynamic test data in an attempt to verify the existence of an empirical 
law to predict structural damping in large aerospace structures. Based 
on a number of previous studies, it was determined that the quantity D, 
energy dissipated/cycle, was the most appropriate measure of damping. 

A number of different hypotheses were checked. It was found that for 
bending vibrations, the formula 
D = 0. 286 T 0 * 746 

where T is the peak kinetic energy and units are in Newton-meters 
predicted damping to within + 2 db levels for 60% of the data examined. 
Since more than half of the data pertained to one vehicle, another formula 
D = 0. 153 T 0 - 893 

was developed to fit the rest of the data within a + 2 db range for 71% of 
the data points. 

Formulas for torsional and longitudinal vibrations were also 
found. It was found for these cases too that the best law, in the sense 
of a least -squares fit, related dissipated energy to kinetic energy as an 
independent parameter. The formulas found, however, differ from the 
formula for dissipated energy in bending in that the exponents are greater 


30 



than unity. Therefore, at higher energy levels, damping will increase, 
in contrast with bending. A relatively small number of data points were 
available for torsional and longitudinal tests, so that the validity of con- 
clusions established for these cases is indeterminate. 

There is a considerable amount of scatter in the D vs. T data. 

The error associated with the various dynamic parameters was estimated 
and an analysis made to determine the probable error in D and T. It was 
found that a good part of the scatter, but still not all, could be due to 
measurement, data acquisition, and numerical error. 

Because of the enormity of the testing program required to come 
up with just one data point, it is an unfortunate fact that not nearly enough 
data are available for a reliable statistical analysis to determine a trust- 
worthy damping law. This situation is complicated by the use of different 
test procedures, equipment, and numerical methods, all of different de- 
grees of accuracy, in different systems of units, and with different de- 
grees of care in documenting. The present study, then, cannot do more 
than establish a quantitative relationship of encouraging, but not definitive, 
reliability. 

It is recommended that the data from the recently performed 
Skylab modal survey be analyzed in the manner of this study to add to the 
statistical base for determination of an empirical damping law. 
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As an aid to future researchers of structural damping in large 
structures, all the raw data used for the analyses in this report are given 
in Appendix A, all converted to the Standard International system of units. 
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Table 1 

Dynamic Tests for Damping 


— 

1 

Vehicle Description 

Type of 
Test 

No. Data 
Points 

Ref. 

1 

Saturn I 

Bending 

196 

1 

2 

Saturn V, Configuration I* 

Bending 

19 

6 

3 

ii n 

Torsion 

10 

7 

4 

it " 

Longitudinal 

11 

8 

5 

Saturn V, Configuration II** 

Bending 

15 

9 

6 

II M 

Torsion 

6 

10 ■ 

7 

II *1 

Longitudinal 

21 

11 


Saturn V, Configuration II, 




8 

(MSFC C. O. 201)*** 

Bending 

14 

12 

9 

ii n 

Torsion 

3 

13 

10 

n ii 

Longitudinal 

7 

14 


S-IVB* 

Bending 

26 

15 

12 

II II 

Torsion 

6 

15 

13 

S-IVB* SAD-202 

Bending 

20 

16 

14 

u it 

Tor sion 

9 

16 

15 

Boeing 747 

Bending 

16 

5 


* Configuration I consists of the entire Apollo Saturn V vehicle 

** Configuration II consists of the S-II stage, S-IVB Stage, instrument 
unit, and the Apollo spacecraft 

*** a modified configuration II for a stronger torsional support at the 
command module / service module interface. 

$ These configurations consist of the S-IVB Stage, instrument unit, 
and Apollo spacecraft 
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Table 2 


Damping Laws for Bending 

The results given include the data 
from test nos. 1, 2, 5, 8, 11, 13 
and 15 (330 data points) 


Damping Law 

Percentage of data points within Hh 1, +_ 2, 
+ 3, + 4 db band 

+ 1 db 

+ 2 db 

f 3 db 

+ 4 db 

” D =.286 r' 4 ° 

32% 

60% 

80% 

88% 

*D = 261 X* 371 

21 

39 

50 

63 

D = 694 Co" 1 ’ 04 

18 

31 

41 

51 

D = 2. 82 m"° 091 

11 

18 

24 

31 

. 593 1. 020 

* D = 4. 31 m X 

33 

55 

73 

82 

, .020 -2. 15 

D = 1485 m oo 

17 

30 

41 

51 

* D = 270 X- 685 M - 202 . 

21 

37 

51 

63 

* D =. 785 r 678 (X/L)- 1 0828 . 

33 

1 

60 

79 

87 

. 675. 1.04 „1, 369 
* D = .671 m to X 

35 

62 

79 

87 

, . . 620 . 646 / v - / T x 1 • 244 

*D= 261 m to (X/L) 

32 

56 

73 

84 

D = . 596 + .095 T - .000015 T 2 

11 

28 

60 

66 


❖ Does not include test 15 (X is not available from this test) 

Dimensions are in SI units as follows: 

D, T (Newton-meters), X( meters), oo (radians /second), m (kilograms) 


> 
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Table 3 


Damping Laws for Bending 
(without Saturn I data) 

The results given include the data from 
tests 2, 5, 8, 11, 13, 15 (174 points) 


Damping Law 

Percentage of data points within 
+ 1-4 db band 


+ 1 db 

+ 2 db 

+ 3 db 

+ 4 db 

D = .153 T* 893 

45% 

71% 

89% 

95% 

\ 

*D = 63. 4 X- 190 

22 

38 

53 

66 

D = 225 647 

16 

32 

50 

60 

D = 16. 2 m“ * 110 

15 

25 

39 

52 

*D = 3. 41 m* 626 X X - 015 

37 

63 

79 

88 

D = 460 m* 012 - 1 * 42 

15 

30 

49 

64 

*D = 294 X* 082 -1 * 103 

15 

35 

51 

63 

*D = . 132 T ,882 (X/L)" > 014 

47 

75 

92 

96 

*D = . 080 m* 881 1,731 X 1,741 

47 

75 

92 

96 

*D= 300 m' 835 • 807 (X/L) 1 - 529 

34 

58 

82 

90 


*Does not include test 15 (X is unavailable from this test). 


Dimensions are in SI units as follows: D, T (Newton-meters), X (meters), 
(radians/second), m (kilograms). 
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Table 4 

Damping Laws for Torsion 

The results given include the data from 
test nos. 3, 6, 9, 12 and 14 (26 data points) 



Percentage of Data Points within jf 1, +2, 3, _+ 4 db 

band of basic law 

Damping Law 


+ 1 db 

+ 2 db 

+ 3 db 

+ 4 db 

D..101 T 1 ' 279 


33% 

66% 

85% 

96% 

D = 870 X' 378 


4 

22 

48 

67 

D = 13.4 «“* 415 


4 

15 

30 

33 

D = 95.7 m"* 405 


19 

26 

33 

52 

. 946 1. 977 

D = 245 m X 


22 

37 

59 

81 

_ . 0/1 -.476 1.094 

D = 4. 94 m a) 


22 

33 

37 

56 

d. .587 x- 998 „ 2 - 431 


22 

52 

56 

67 

D = .934 T 1-U (X/L)' 150 


29 

74 

93 

100 

_ - Q _ 1.01 2. 69 „2. 33 

D = . 095 m a) X 

37 

81 

93 

100 

^ . 588 2. 573. . ,1. 456 

D = 4. 77 m u (X/L) 

41 

78 

96 

100 

D = -.524 + . 333 T - .0013 

T 2 

37 

59 

63 

70 


Dimensions are in SI units as follows: D, T (Newton-meters), 

2 

X (radians), co (radians /second), m(kilogram-meter s ) 
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Table 5 


Damping Laws for Longitudinal Motion 

The results given include the data from test nos. 
4, 7 and 10 (39 data points) 



Percentage of Data Points within _+ 1, 
+ 4 db band of basic law 

+ 2, +3, 

Damping Law 


+ 1 db 

+ 2 db 

+ 3 db 

+ 4 db 

D = .057 T 1 ' 104 


41% 

67% 

90% 

100% 

D = 269 X* 192 


13 

23 

44 

51 

D = 235, 930 to" 1 ’ 38 ° 


15 

38 

49 

56 

D = .993 m‘* 191 


13 

26 

38 

49 

D = 25. 1 m* 847 X 


38 

59 

72 

77 

D = 42, 108 m #159 co'* 262 


18 

31 

49 

54 

-D= 590, 740 X- 158 Co" 2 - 433 


18 

36 

51 

56 

D = .0084 T 1,158 (X/L)”° 


49 

74 

82 

100 

1. 259 3. 315 2. 405 

D = . 0002 m to X 

54 

77 

92 

100 

„ co ^ 1. 164 2. 379. /T ,2. 145 

D = 58. 37 m to (X/L) 

26 

59 

79 

82 

( 


Dimensions are: D, T (Newton-meters), X (meters), co( radians /second), 
m (kilograms) 
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Table 6 


Comparison of Dissipated Energy Values 


Test No. 

No. Data Pts. 

Average percent difference for 
D calculated two ways* 



Bending . 


3 

16 

73. 6 

4 

18 

91.1 

5 

19 

13. 1 

6 

14 

1.0 

19 

15 

12. 0 


t 

Torsion 


9 

3 

65. 1 

10 

6 

37. 3 

12 

5 

50. 3 

13 

3 

24.8 

17 

10 

42.8 


Longitudinal 


15 

21 

\ 

\ 

315. 6 

18 

11 

\ 

J 

s' 

9.0 

21 

7 

247.8 


*The average percent difference was calculated for each test by 
taking the average value of the quantity 100 x ID } - D 2 I/D 2 
where D = it FX e and D = 47r f T. For each case, a few 
data points showing untypically large error were excluded 
from the average. 
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Figure 2, Dissipated Energy/Cycle vs. Maximum 
Kinetic, ^nergy for Bending Tests 
(Exclud_ .g Saturn I Data) 
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Figure 3. Dissipated Energy/Cyc^’ vs. Maximum Kinetic Energy for Torsional Te^ 1 
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Figure 4* Dissipated Energy /Cycle vs. Maximum Kinetic Energy for Longitudinal Tests 


APPENDIX A 


DYNAMIC TEST DATA USED FOR DETERMINATION 
OF EMPIRICAL DAMPING LAW 







test 

NO. 1 SATURN 

I DYNAMIC 

TEST VEHICLE SAD^6 PITCH 

+ YAW 







FRECUWY 

FORCE 

XN 

XE 

gene mas 

s 

^PzETA 

* 

DC 

= PI*XE>F] 


T 

D C = 

4*PI*|P 

^T3 


[HZ] 

CNl 

CKl 

CM] 

[KG-SEC2/M] 




CN-M] 


CN-M] 


CN-V] 


* 

X 

3 

•57o 

2284. 9A9 

1 »6800E*-02 

1 .7741E-03 

8* 1984E 

0? 

. OOOOE 

00 

1 

♦2735E 01 

5 

•8212E 

01 

•OOOOE 

,00 

2 

8 

• 860 

2088.816 

6.3000E-04 

1 .8648E-04 

2 • 6478E 

03 

•0000E 

00 

1 

• 2237E 00 

1 

• 6284E 

00 

•OOOOE 

00 

3 

14 

• 44q 

1117.958 

1.7100E-04 

4.859SE-05 

9.12C2E 

0? 

♦OOOOE 

00 

1 

•7068E-01 

1 

•0976E 

-01 

-OOOOE 

00 

4 

1 

•600 

2098.623 

3.4650E-03 

5.53C1E-04 

3 « 2470E 

04 

•OOOOE 

00 

3 

* 64 60E 00 

1 

• 9700E 

01 

•OOOOE 

00 

5 

2 

• 450 

2083.816 

6.5500E-04 

1 .9034E-04 

7.3795E 

04 

•OOOOE 

00 

1 

• 249 1 E 00 

3 

• 75 1 2E 

00 

•OOOOE 

00 

6 

3 

•220 

2128.043 

5.7180E-03 

2.8018E-04 

1-8927E 

03 

•OOOOE 

00 

1 

* 3 7 3 1 E 00 

1 

•2665E 

01 

•OOOOE 

CO 

7 

3 

•300 

2235.916 

1 • 1090E-03 

2.76l4r- 0 5 

6 • 2076E 

03 

•OOOOE 

00 

1 

•9397E-Q1 

2 

• 1 76 1 E 

00 

•OOOOE 

00 

8 

4 

•29o 

2235.916 

5.930CE-04 

2.C126E-Cf4 

1-5063E 

04 

•OOOOE 

00 

1 

» 4 1 3 7E 00 

1 

. 9243E 

00 

•OOOOE 

00 

9 

5 

•020 

2118.236 

3*52SoE^o3 

1 .9369E-04 

2.392SE 

03 

•OOOOE 

00 

1 

•2883E 00 

1 

• 4 S i 5 E 

01 

• 0090- 

00 

10 

5 

•800 

2118.236 

3 • 9200E-05 

1 .2948E-05 

1 . 8S39E 

04 

•OOOOE 

00 

8 

■ 6 1 63E-02 

1 

•9222E 

-02 

•COOOE 

00 

11 

7 

• 2 1 1; 

2118.236 

2 • 01 CCE-05 

6.6792E-06 

2.5890E 

04 

•OOOOE 

00 

4 

.444SE-02 

1 

• 0733E 

-02 

•OOOOE 

00 

12 

7 

• n 00 

2265.336 

1 «960CE-04 

7.1266E-05 

2 • 2?20E 

04 

•OOOOE 

00 

5 

•0718E-01 

9 

• 4 759E 

-01 

• OOOOE 

00 

13 

10 

• 49o 

2177.076 

3«9000E-05 

1 .1376E-05 

2-8204E 

04 

•OOOOE 

00 

7 

.7808E-02 

9 

•3179E 

-02 

•CO "SOE 

00 

1* 

13 

• S7o 

2128-043 

7 *0000E-05 

2.5263E-05 

3-46o8E 

04 

•OOOOE 

00 

1 

•68S9E-01 

6 

• 4395E 

-01 

• COCOE 

00 

15 

1 

• 720 

2059.396 

1.9072E-02 

1 .2633E-03 

3.2950E 

03 

•OOOOE 

00 

8 

•2055E 00 

6 

• 999oE 

01 

•CO''OE 

00 

16 

1 

•600 

2137.850, 

3»30 6 OE-o3 

5.431 8 E”04 

3 • 2382E 

04 

•OOOOE 

00 

3 

• 64 8 1 E 00 

1 

• 7S3->E 

01 

•OOOOE 

00 

17 

1 

•73o 

2128.043 

1 . 8 1 4 2E-Q2 

1 .2845E-03 

3.3048E 

03 

•OOOOE 

00 

8 

« 587 1 E 00 

6 

• 426QE 

01 

•OOOOE 

00 

IS 

2 

• 45 0 

2118.236 

6-2500E-04 

1 .9562 E -o4 

6.80SSe 

04 

•OOOOE 

00 

1 

•30l?E 00 

3 

• 1 513E 

00 

•OOOOE 

00 

19 

3 

.230 

2128.043 

5 . 7 1 70E-03 

2.8128E-04 

1 . 951 5e 

03 

.OOOOE 

00 

1 

,8805e 00 

1 

, 3 1 35e 

01 

• OOOOE 

CO 

20 

3 

• 810 

2157.263 

8.9400E-04 

3.0575E-05 

9-2034E 

03 

•OOOOE 

00 

2 

•0723E-01 

2 

• 10S8E 

00 

•OOOOE 

00 

21 

4 

• 370 

2157.463 

3 • 3900E-Q4 

1 .3723E-04 

1 .8642E 

04 

•OOOOE 

00 

9 

•301 1E-C1 

8 

•0760E 

-01 

•OOOOE 

00 

22 

5 

•0^0 

2137.850 

3.5P10E-03 

2 .0492E-04 

2 .43 HE 

03 

•OOOOE 

00 

1 

•3763E 00 

1 

• 54SLE 

01 

•OOOOE 

00 

23 

5 

*810 

2128.043 

4.6300E-Q4 

1 .4955E-04 

1 • 61 q3E 

04 

•OOOOE 

00 

9 

• 9980E-0 1 

2 

OOOOE 

oc 

•CO^OE 
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3336.165 

1 

•5734E-Q2 

6 

•4010E-04 

2* 

9925E 

03 

9 

•0000E-03 

7* 

461 IE 

00 

6 * 59 / 1 E 

01 

7*461 jE 

00 

10 

2*976 

16000*847 . 

9 

* 44 3q2-o3 

6 

•0532E-04 

9* 

1000E 

03 

1 

♦2000E-02 

2* 

1 392E 

01 

1.4186E 

02 

2 • 1 392E 

01 

11 

1 *S£ u 

5386.794 

•3 

•9926E-02 

2 

•0643E-03 

4* 

4450E 

03 

7 

-0000E-03 

2* 

9942E 

01 

3 « 4 038E 

02 

2*9942E 

01 

12 

4 * ^ 6 Q 

2 1 6o5 • CO^ 

5 

• 1 98 1 E -q3 

8 

•1S83E-04 

1 ■ 

l^oOE 

O'* 

9 

•0000E-03 

1* 

4278E 

01 

1«2625E 

02 

1 • 4278E 

01 

13 

1*112 

4359.256 

5 

•78SQE-02 

4 

.7642E-03 

1* 

1 357E 

O'* 

6 

•0000E-03 

7* 

0022E 

01 

9.2S69E 

0? 

7*0Q^2E 

01 

1* 

1*810 

334Q.513 

2 

*0'^3E-c2 

8 

♦313SE-04 

3* 

78 0 0E 

03 

7 

•5000E-03 

9. 

627 SE 

00 

1 • C215E 

02 

9 • 62 7 8E 

00 

15 

1*821 

4159.085 

2 

• 4 4 44E - q2 

1 

. 0296E-03 

3* 

64poE 

03 

8 

«5000E-03 

1* 

5207E 

01 

1 -4237E 

02 

1 • 52p7E 

01 

1^ 

2*578 

6690*123' 

9 

• 8936E - 03 

3 

•09l8E-04 

3* 

5350E 

03 

9 

•0000E-03 

5* 

1339E 

00 

4.5394E 

01 

5 • 1 339E 

00 

17 

3 9 45q 

21702.865 

4 

•3377E-03 

3 

.7675E-04 

4 • 

1475E 

O 4 * 

1 

♦2000E-02 

2 * 

7649E 

01 

1 .8335E 

02 

2 • 764 9E 

01 

18 

1*258 

3385.095 

5 

•2200E-02 

4 

.3360E-03 

6 • 

7200E 

0^ 

6 

•0000E-03 

4 » 

3127E 

01 

5-7200E 

02 

4 • 3i 37 F 

Cl 

19 

2*138 

435Q.359 

1 

•8S67E-02 

7 

•7596E-04 

2* 

8875E 

03 

8 

♦0000E-03 

9* 

323SE 

00 

9 • 2 7 4 6 E 

01 

9 • 3H.38E 

00 


•*! 

D 



TEST N0- 3 SATURN V OTV C0NEIG- I TORS 1 0NAL (1967) 



FRFOulMr 

T0RQUE 

XN 

XE 

genl mass 

^2ETA 


DC«PI*XE*F) 

T 

DC- 

4*Pl *^^Ut3 


CHZ3 

[N-M] 

craD3 

CRAD3 

CKG-M23 




CN-MJ 

CN-MD 


cn-mj 

1 

5 • 66o 

24756-178 

6-0920E-05 

5-4780E-06 

1-0620E 

06 

1-0000E-02 

i 

3-1 320E-0 1 

2 • 4 924 E 

00 

3- 1320 e "-01 

2 

6.162 

25650-979 

8-469QE-05 

5 -9030E-06 

3-38942 

05 

S-OOOOE-03 


1.8317E-01 

1 -8220E 

00 

1 *83172-01 

3 

8*117 

23115-70S 

3 - 1 920E -05 

4.9104E-06 

2-8245E 

06 

1 -6000E-02 


7.5252E-01 

3-7427E 

00 

7-5252E-01 

4 

8-339 

20580-437 

1 -4860E-04 

8 • 1 972E -06 

4-3836E 

05 

1 -0000E-02 


1-S759E 00 

1 -4928E 

01 

1-8759E 00 

5 

11-326 

27142.315 

6-7q30E-05 

1 -9219E-05 

3 • 05q5E 

05 

2 -OOOOE -02 


8.7222E-01 

3-4705E 

00 

8-7222E-01 

6 

5-691 

22708-980 

6-2710E-05 

6-0300E-06 

1 • 028 1 E 

06 

9-0000E-03 


2-9233E-01 

2-5848E 

oc 

2-9233E-01 

7 

6-194 

22708-980 

6 • 4320E-05 

6-0300E-06 

7-6826E 

05 

9 -OOOOE -03 


2-7222E-01 

2-4070E 

00 

2 • 7222E- 0 1 

8 

8-113 

21918-120 

3-2980E-05 

7 -7600E-06 

3-5024r 

06 

1-4000E-02 


8-7075E-01 

4 - 9494E 

00 

8-7075E-01 

9 

8-895 

23160-500 

1 -4430E-04 

2.8700E-05 

4-5192E 

05 

1 -1000E-02 


2-03i5e 00 

1-4697F 

01 

2*03l5E 00 

1C 

11.423 

26389.240 

6-6400E-05 

1 .1900E-05 

3.8A13E 

05 

1.8000E-02 


9 - 8930E-01 

4-3737E 

00 

9.8930E-01 


°1 

IS** 



TEST IM0. 4 SATURN V 9TV C6NEIG. I LSNG I TUD I NA.L (1967) 



FReQU^^T 

F6RCE 

XN 

XE 

genl mass 

^Pzeta 


DC=PI*XE 

*F1 

T 

D C = 4*P I T] 


CH.2] 

CN3 

CM] 

CM] 

CKG-SEC2/M] 



r~i 

r 

i 

? 

Ui 


CN-M] 


Cn-*0 


1 

3*752 

765 0 9.384 

1.8547E-03 

9*275lE-0' t 

2«2400E 

0* 

8.0000E-03 

f 

2.1526E 

02 

2* 1 4 1 2E 

03 

2 ♦ 1 5 ? 6 E 

02 

2 

4.463 

78955.905 

4 • 2303E -03 

1 -3114E-03 

3 • 92C0E 

05 

9-OOOOE-03 


3.1193E 

02 

2 . 758 1 E 

03 

3-1193E 

02 

3 

6 • 506 

84516.1 80 

4*1 736E -04 

8.7641 E”05 

9.5375E 

05 

1 .3000E-02 


2.2676E 

01 

1 .3881E 

02 

2-2676E 

01 

4 

7.568 

77843.850 

6-8362E-0f> 

1 -0937E-04 

2 * 9750E 

05 

1 -4000E-02 


2.7654E 

01 

1 -5719E 

02 

2 • 76" 4E 

01 

5 

5.114 

78066.261 

3.Q855E-C3 

9.8732E-04 

2-6250E 

03 

1 .7000E-0? 


2.7560E 

02 

1 .2901E 

03 

2 • 75 fc oE 

02 

6 

6-492 

69837.054 

1.8720E-03 

2.2467E-04 

2-9750E 

05 

6 • 0000E-03 


6 * 5393E 

01 

8 . 673 1 E 

02 

6 »5353E 

01 

7 

8.276 

68280*177 

8.9686E-05 

3.6771E-04 

4.3750E 


1 -3000E-02 


7.7723E 

01 

4.7577E 

02 

7 • 77?3E 

01 

8 

9*376 

69837.054 

9-1699E-05 

7 •0827e*0‘5 

7 . 0000E 

06 

l*0000E-02 


1 .2S35E 

01 

1 • 05 1 4 E 

0? 

1 • 2835E 

01 

9 

6» 176 

56714. 8q5 

3-Q540E-03 

5 •4969E-0'* 

1*1 637 E 

05 

9*0000E-03 


9.2426E 

01 

8* l ’22E 

02 

9 • 24 ?6E 

01 

10 

7.976 

53378.640 

3.95S9E-0 1 * 

2.4544E-Q3 

1 . 0&75E 

07 

1 -4000E-02 


3.6962E 

02 

2* 1010E 

03 

3*696?E 

02 

11 

9.406 

64054.368 

1 -6241E-Q3 

3.2482E-0'* 

6.8250E 

0* 

1 "0000E-02 


3 « 95q8E 

01 

3. 1439E 

02 

3 • 9 5 n 8 E 

01 


f\) 



TEST NS < 5 SATURN V DTV CSNF 



FReQuWP' 

FSRCE 

XN 

XE 

genl mass 


[HZ3 

Cn3 

CM] 

CM] 

CKG-SEC2/M] 

1 

1 • 667 

6325.369 

4.9193E-02 

1 .4745E-03 

3.3425E o3 

2 

2*701 

6707.916 

9.6C19E.-03 

1 .9204E-04 

3-7450E 03 

3 

5*979 

9194.47] 

7 » 0235E-C'* 

1 •2639E-0'» 

1-43S5E 05 

4 

2-224 

3927.778 

6 * 20®4E - o3 

8.1153E-0'* 

3 • 5S75E o'* 

5 

5*953 

9425.778 

4.Q995E-04 

1 - 106 &E- 0 A 

4*2350E 03 

6 

8*555 

8660*684 

2-1 146E-Q5 

8*855lE-0^ 

2 ♦ 8000E 03 

7 

2*229 

4870*801 

7.6645E-03 

9 . 1974 E-0^ 

3.9025E O'* 

8 

5*897 

17321.369 

7 • 1487E-04 

7.14S7E-05 

3.2900E 05 

9 

8*690 

5791.582 

l«0?36E-05 

5.3536E-04 

4 . 7775E 08 

10 

2*264 

2993.652 

8.6757E-03 

2.1686E-03 

3.5175E o'* 

11 

5*916 

14950.467 

■ 4.2223E-0A 

1 .0558E-04 

5.7400E 05 

1? 

1*66 1 

5489.103 

3 . 42 10E-02 

6.8430E-0'* 

3.8150E 03 

13 

2*701 

127o4.1l6 

1.4509E-02 

4 *35i7e;-o4 

4 . 3575£ 0 3 

1* 

5*733 

7241 .7q2 

2 • 3669E-0 1 * 

1 «0l75E-0 i * 

5.5125E 05 

15 

6*303 

8416.Q32 

2.2004E-0'* 

1 • 2S95E-0'* 

7.7875E 05 


Ga^Il PITCH * YAW FREO. 

SWEEP (1967! 





• 

^ETA 

D C = PI»XE*F] 


T 

DC 3 

4«PI 

g 

CN-M] 



CN-M2 



5«2000E-03 

2.8993E 

01 

4 

* 4368E 

02 

2*89932 01 

7.4000E-03 

4 » 6 236E 

00 

4 

• 9721 E 

01 

4*62362 00 

6*3000E-03 

3 * 9642E 

00 

5 

•0073E 

01 

3 * 964 ?E 00 

6«3000E-o3 

1 .q588E 

01 

1 

• 35 0 1 E 

02 

1*q 6S8E: oi 

6«3000E-03 

3* 94 l5E 

00 

4 

» 9787E 

01 

3 * 94 i5E 00 

1.0300E-02 

2.3413E 

01 

1 

• 808SE 

02 

2 • 34 i 3E 01 

6.3000E-03 

1 .7SO0E 

01 

2 

•2483E 

02 

1 *78^0- 01 

7.3000E-03 

1*0587E 

01 

1 

• 1 54 1 E 

02 

1*C*587E ot 

9 • 5000E - 03 

9 . 437 1 E 

00 

7 

• 9o5 IE 

01 

9 * 4 3 7 l E oo 

5.8000E-03 

1.9524E 

01 

2 

• 6737E 

02 

1*952 4 L ol 

9-0000E-03 

7.9955E 

00 

7 

•0695E 

01 

7*99=52 00 

6.0000E-03 

1.8332E 

01 

2 

.'♦S^E 

02 

1*83322 01 

9 • 3000E -03 

1 .54382 

01 

1 

• 3210E 

02 

1*5 4 3 8 E 01 

9.6000E-03 

2.4162E 

00 

2 

•0029E 

01 

2*41 62E 00 

9.0000E-03 

3*3‘*41E 

00 

2 

• 9568E 

01 

3 • 344 1 E 00 









test 

K8. 5 SATURN 

V 

DTV C6NEIG. II PITCH + 

YAW FORCE LEVEL ( 1967) 





i 


frecu/Bt 

FORCE 


XN 

XE 

6ENIL MASS 

WEETA 


DC 

=PI*XE 

*F] 


T 

DC* 

4* p I * zl^ 

wn 


CHZ] 

C\'3 


CM] 

CM] 

CKG-SECc 

‘/Ml 




CN-M] 



CN-M] 


CN'-M] 


1 

1-665 

7210*565 

5 

•5195E-02 

1.5978E-03 

3* 

34258 

03 

•OOOOE 

00 

4 

3 

.6193E 

01 

5 

• 5722E 

02 

•00009 

00 

2 

1 • 67 1 

5689.273 

4 

•7334E-02 

1 .3995E-03 

3- 

3425E 

03 

•OOOOE 

00 

2 

• 501 “E 

01 

4 

* 1 276E 

02 

•oonoE 

CO 

3 

1*676 

3829.917 

3 

•7385E-02 

1 *1193e-03 

3* 

34258 

03 

•OOOOE 

00 

1 

*3467 E 

01 

2 

•59o3E 

02 

•CCOOE 

00 

4 

2*710 

6965.913 

9 

8521E-03 

2.5743e-0“ 

3* 

74508 

03 

•OOOOE 

00 

5 

.63358 

00 

5 

• 26968 

01 

•OOOOE 

00 

5 

2*720 

5564.723 

7 

•9806E-03 

2»1002E-0“ 

3* 

74508 

03 

•OOOOE 

00 

3 

.67158 

00 

3 

•48338 

01 

• COOOE 

00 

6 

2*722 

3696.471 

5 

•9621E-03 

1.61598-0“ 

3. 

7450E 

03 

•OOOOE 

00 

1 

*87658 

00 

1 

• 947QE 

01 

•OOOOE 

00 

7 

2*224 

3491.853 

5 

.2434e-o3 

7 .09638-0“ 

3. 

53758 

04 

•OOOOE 

00 

7 

.78478 

00 

9 

.62988 

01 

•OOOOE 

DO 

8 

2*223 

2388*694 

3 

•6171E-03 

“ .966SE-04 

3. 

5S75E 

04 

♦OOOOE 

00 

3 

• 7-72E 

00 

4 

•5785E 

01 

•OOOOE 

QO 

9 

2*225 

1797.081 

2 

•7366E-03 

3 * 7097E *04 

3. 

5375E 

04 

* OGOOE 

00 

2 

• 09““E 

00 

2 

• 6? J5E 

01 

•OOOOE 

00 

1C 

5*955 

5395.691 . 

2 

•5661E-04 

7.0439E-05 

4 . 

2350E 

05 

•OOOOE 

00 

1 

. 1940E 

00 

1 

•9520E 

01 

•COOOE 

00 

11 

5*9F.t 

3647.540 

• 6 

•9479E-05 

4*961 1 E-05 

4. 

2350E 

05 

•OOOOE 

00 

5 

• 6 850E 

-01 

1 

• 4315E 

00 

•OOOOE 

00 

12 

5*954 

2771.241 

4 

•7147E-05 

3.5753E-05 

4. 

2350E 

05 

•OOOOE 

00 

3 

.11278 

-01 

6 

• 587-.E- 

-01 

•C0009 

DO 

13 

8*765 

9999.599 

1 

* 4247E -05 

4.4S48E-04 

2* 

8000E 

O'? 

•OOOOE 

00 

1 

.40398 

01 

8 

. 6 1 8 1 E 

01 

• Cl oc9 

00 

14 

8*766 

7139.393 

1 

•5306E-05 

4.1248E-0'* 

2. 

8000E 

os 

•OOOOE 

00 

9 

• 25158 

00 

9 

• 3495E 

01 

•OOOOE 

00 

15 

8*764 

5422*380 

9 

• 1464E-Q6 

3 .59948-0“ 

2* 

8000E 

os 

•OOOOE 

00 

6 

•1316E 

00 

3 

• 55 1 3E 

01 

•OOOOE 

00 

16 

2*222 

5449*069 

8 

•6902E-03 

1 • 3 1 28E *03 

3. 

9q25e 

04 

•OOOOE 

00 

2 

•2“738 

01 

2 

• 8722E 

02 

•OOOOE 

00 

17 

2*223 

3629.748 

5 

•6194E-03 

8 .4QS7e-0“ 

3. 

9C25E 

04 

•OOOOE 

00 

9 

• 58858 

00 

1 

•2C21E 

02 

♦OOOOE 

00 

18 

2*224 

2620*002 

4 

.2596E-03 

6. 35988-0** 

3. 

9025E 

04 

•OOOOE 

00 

5 

.23“8£ 

00 

6 

.9132E 

01 

•OOOOE 

00 

19 

5*876 

17023.338 

7 

•2412E-0“ 

1 .5143E-0“ 

3. 

2900E 

05 

•OOOOE 

00 

8 

>09878 

00 

1 

• 1757E 

02 

"OOOOE 

00 

20 

5*883 

12566.221 

5 

•7989E-C“ 

1-1578E-Q4 

3. 

2900E 

05 

•OOOOE 

00 

4 

•57o6E 

00 

7 

• 55828 

01 

•OOCOE 

00 

21 

5*898 

8656.236 

4 

•1411E-0'* 

8.6573E-05 

3 • 

2900E 

09 

•OOOOE 

00 

p 

.35438 

00 

3 

•874 0 E 

01 

*00009 

00 

22 

8*620 

10 4 08*F.35 

1 

*8227E-o5 

1 .2329E-03 

4. 

77758 

0 ? * 

•OOOOE 

00 

4 

•19538 

01 

2 

* 328qE 

02 

•C0009 

00 

23 

8*621 

7l7 C .53l 

1 

. 3954E-0 1 ’ 

7 .75848*0“ 

4 • 

77758 

os 

•OOOOE 

00 

1 

.74778 

01 

1 

•36488 

02 

•00009 

00 

2 4 

2*268 

3629.748 

9 

*4575E-o3 

2*3526 E - 0 3 

3* 

51758 

04 

•OOOOE 

00 

2 

.68278 

01 

3 

•1945E 

02 

* OOCOE 

00 

25 

2*27o 

' 1712.565 

4 

• 4 ].72E-q3 

1 • o7S3e-o3 

3. 

51758 

O 4 

•OOOOE 

00 

5 

.3015E 

00 

6 

• 9S 0 9E 

01 

•OOhqE 

00 

26 

5. 905 

1479q.331 

4 

• 7674E-0'* 

8 *9o38e-o5 

5. 

74Q0E 

05 

•OOOOE 

00 

4 

.13728 

00 

9 

. 1013E 

01 

•OOCOE 

00 

27 

5*946 

1 1285.134 

3 

•4224E-o“ 

6,1675e-o5 

5. 

7400E 


•OOOOE 

00 

2 

.18668 

00 

4 

•69208 

01 

• co^qe 

00 

28 

5*943 

7788.833 

2 

•53698-0“ 

3*20l9E-05 

5. 

74oOE 

05 

-OOOOE 

00 

7 

.83498' 

-01 

2 

• 5755E 

01 

•CO^flE 

00 

29 

1*658 

5217.762 

3 

•9325E-03 

1 .1660E-03 

3. 

8150E 

O 3 

•OOOOE 

00 

1 

• 91128 

01 

3 

• 2C13E 

02 

• 0 0 - 0 E 

CO 

30 

1*653 

4328. 1 13 

3 

•5609E-02 

1 *o375e-0 3 

3 * 

8150E 

03 

•OOOOE 

00 

1 

• 410?E 

01 

2 

•6249E 

02 

•OOOOE 

00 

31 

1 • 6? 8 

3487. 4 0 4 

2 

•9q“6E-o2 

8*14728-0“ 

3. 

8150E 

03 

•OOOOE 

00 

8 

•9260E 

00 

1 

•7“658 

02 

•OOCOE 

00 

32 

2* 7q1 

12761.943 

1 

.4946E-02 

29538-0“ 

4. 

3575E 

0 3 

•OOOOE 

00 

1 

.72338 

01 

1 

» “ 0 1 8 E 

02 

, ' *00009 

00 

33 

2* 7q1 ■ 

93q5 * 676 

1 

*2181E-c3 

3 . 43328 - 0 “ 

4. 

3575E 

03 

•OOOOE 

00 

1 

•01838 

01 

9 

• 3 1 0 IE 

01 

L. . U — 

00 

34 

2*701 

6338.713 

6 

*9S7cE-03 

1*92148-0“ 

4. 

35758 

0’3 

•OOOOE 

00 

3 

. 8262E 

00 

3 

•0633E 

01 

•OOOOE 

DO 

35 

5*75 0 

6752. 39S 

2 

•“0 65 E-0“ 

8.80“7e-05 

5. 

51258 

05 

•OOOOE 

00 

1 

.86788 

00 

2 

•08348 

01 

•OOOOE 

00 

36 

5.75 0 

5137.694 

1 

• 32778-0“ 

7 *0550E-05 

5. 

51258 

05 

•OOOOE 

00 

1 

.13878 

00 

6 

* 34 1 8E 

00 

•COOOE 

00 

37 

5*75o 

2219.662 

1 

•l735E-o“ 

2.9S438-05 

5. 

5125E 

05 

•OOOOE 

00 

2 

>0 6 OlE 

-01 

4 

. 954q£ 

00 

•OOOOE 

00 

38 

6*29q 

7993.451 

2 

•“37qE-0“ 

1 * l73“8-0“ 

7. 

78758 

05 

•OOOOE 

00 

2 

.94668 

00 

3 

•6120E 

01 

•OOOOE 

00 

39 

6 • 29q 

6467.712 

1 

•8q75E-o“ 

9.67688-05 

7. 

7875E 

05 

•OOOOE 

00 

1 

. 9662E 

00 

1 

• 9 8 7 0 E 

01 

•00009 

CO 

40 

6*29q 

4412.634 

3 

•992“E-0“ 

6.66438-05 

7. 

7875E 

05 

•OOOOE 

00 

9 

.23858 

-01 

9 

. 6937E 

01 

♦OOOOE 

00 
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TEST 

NO. 6 SATURN 

V DTV ce^Fl 

:g. i i torsional freg. 

SWEEP (1967) 




TORQUE 

CN-M] 

XN 

CRAD] 

XE 

[RAD] 

genu mass 

CKG-M23 

^ZETA 

D [ = PI*XE*F] 
cn-m] 

T 

CN-M] 

DC- 

4*PI#!B\4T3 

cn-mi 

5-605 

17715.264 

5* 1 400E*05 

6.5856E-06 

1 • 5252E 06 

4 

5» 1000E-03 

3* 6652E-01 

2 * 4989E 

00 

l # 6ot5^-01 

6-290 

6134.814 

3-3200E-05 

3.63S1E^06 

1-1637E 06 

4«60Q0E-03 

7 • 01 1 7E-02 

1 • 00 l 7E 

00 

5*75^55-0? 

9*736 

9885.750 

1 • 1 400E-0^ 

1 .2862E-05 

2.2596E 05 

6.0000E-03 

3,9947e-01 

5 , 4945E 

00 

4.-14P8E-01 

5-590 

17398.920 

3.4600E-05 

8.0234E-06 

1-6269E 06 

6.5000E-03 

4.3856E-01 

1-1757E 

00 

9 • 6033E-02 

6*25$ 

6281-688 

3*3400E-o5 

3.5584E-06 

7.9086E o3 

6-3000E-03 

7.Q223E-02 

6.8136E- 

■01 

5-3942E-0? 

9-644 

9648.492 

9»3300E-05 

1 -1017E-05 

2.1353E qS 

8-0000E-03 

3.3393E-01 

3.4125E 

00 

3-43f6E-(jl 
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TEST 

NO. 6 SATURN 

V DTV c 

ONF I G . 11 TORSIONAL FORCE 

LEVEL (1967) 



- 



FReoJ^Y 

T8RQUE 

XN 

XE 

GENL MASS 

^PzETA 

* 

D [ = P I *XE#F] 


T 

0 [ = 4*P 1 

A*T3 

CHZ5 

CN-M1 

CRA01 

[RADI 


CKG-M2] 




CN-MD 


CN-M] 


CN-m] 


5 * 605 

17686.703 

l*0200E-0‘» 

7.2276E-06 

1 

•5252E 

06 

•OOOOE 

00 

4.0159E-01 

9 

•84075 

00 

•OOOOE 

00 

5-611 

13649.521 

8.5837E-05 

5.3921E-06 

1 

• 5252E 

06 

•OOOOE 

00 

2.3133E-01 

6 

• 9838E 

00 

•OOOOE 

00 

5.620 

8362.734 

6.5058E-05 

3.2544^-06 

1 

-5252E 

06 

.OOOOE 

00 

8.5500E-02 

4 

-0H48E 

00 

•OOOOE 

00 

6-268 

5983.681 

1 .1 181E-04 

4 .0727E-C6 

1 

- 1637E 

06 

•OOOOE 

00 

7 • 6560E-02 

1 

-12S2E 

01 

•OOOOE 

00 

6-275 

4854.231 

9.5650E-05 

3.4714E-06 

1 

• 1637E 

06 

•OOOOE 

00 

5 * 2939E- 03 

8 

•27505 

00 

•COOOE 

00 

6-284 

3532.534 

7 . 4656E-Q5 

2.6367E-06 

1 

• 1637E 

06 

•OOOOE 

00 

2.9261E-03 

5 

• 0556E 

00 

•eooo5 

00 

9-735 

10573.572 

2»5003E-o5 

1 .0147E-05 

2 

-2596E 

05 

•OOOOE 

00 

3.370SE-01 

2 

-64245 

-01 

•0000^ 

00 

9-745 

7665.840 

1.5619E-05 

7 .5641E-0’6 

2 

•2596E 

05 

•OOOOE 

00 

1.8217E-01 

1 

• 0-9345’ 

-01 

•oonoE 

00 

9-775 

5118.C70 

7 . 6360E -q6 

5.2124E-06 

2 

-2596E 

05 

•OOOOE 

00 

8.3318E-02 

2 

* 4850E 

-02 

- OOOO 5 

00 

5 • 51 Q 

17566.548 . 

3.8Q43E-05 

1 .5353E-05 

1 

• 6269E 

06 

•OOOOE 

00 

8.4730E-01 

1 

•41105 

00 

•0000^ 

00 

5*51/ 

13048.749 

•2.6646E-05 

1 »0072e _ 05 

1 

-6269e 

06 

•OOOOE 

00 

4.1288 E -01 

6 

• 9403E 

-01 

•00005 

00 

5.5^0 

§482 . 889 

3.1 883E-05 

7 .34SSE-06 

1 

• 6269E 

06 

•OOOOE 

00 

1 •9318E-01 

4 

*71985 

-01 

•00005 

00 

6*239 

10044.894 

2 • 0609E-04 

9.2344E-06 

7 

-9086E 

05 

•OOOOE 

00 

2.9141E-01 

2 

• 5808E 

01 

•CC 205 

00 

6*237 

7689.871 

1.7417E-0* 

6.3954E-06 

7 

-9086E 

05 

•OOOOE 

00 

1 »5450E-01 

1 

* 842 1 E 

01 

•0000 E 

0^ 

6*241 

5334.843 

1 •2104E-0'* 

3.2487E-06 

7 

•9086E 

05 

•OOOOE 

00 

5.444 8£ -02 

8 

•9079E 

00 

•000C5 

00 

9*635 

3 1919.300 

1.7264E-0* 

2 .6250E-05 

2 

- 1353E 

05 

•OOOOE 

00 

9.8295E-01 

1 

- 1 662E 

01 

♦0000- 

00 

9*636 

9780.554 

1.50*8E-0* 

2.2343E-05 

2 

- 1353E 

05 

•OOOOE 

00 

6 . 8652E-0 1 

8 

-8625E 

00 

•00005 

00 

9*641 

6079.8Q4 

9.4275E-05 

1 .34S5E-05 

2 

• 1353E 

05 

•OOOOE 

00 

2.5756E-01 

3 

•4 8205 

00 

•OOOOE 

00 








TEST 

NO * 7 SATURN 

V 

DTV CONFIG. 

I I LONGITUDINAL FREQ. 

SWEEP (1 

967) 







FRE3 

• 

FfiRCE 


XN 

XE 

GENL MASS 

^ZETA 

DC» 

RI*XE 

*F] 

T 

D t c 4* 

P I * 


CH2 J 

CN3 


CM] 

CM] 

CKG-SECc 

VM] 




CN-M] 


CN-M] 



Cn-*] 


1 

6* 

213 

5529.137 

5 

• 856qE>04 

5.2705E-05 

7- 

6825E 

04 

3* 

8000E-03 

9. 

5858E 

-01 

2*0074E 

01 

9 

•5858E-01 

2 

a • 

2C5 

21769. 5S9 

8 

• 4435E-04 

1.3510E-OA 

9* 

59qoE 

04 

6 • 

5000E-03 

9* 

340 4 E 

00 

1 * 1 435E 

02 

Q 

•34r.4E 

00 

3 

9* 

846 

1 7299 . 12S 

5 

•3811E-05 

l»6143E-04 

1 • 

1S47E 

07 

1* 

2000E-02 

9 

8995E 

00 

6» 5648E 

01 

9 

•■8995E 

00 

4 

11* 

217 

16867.650 

3 

•297qE-04 

6.5941 E -05 

1* 

3177E 

05 

7. 

6000E-03 

3 

3977E 

00 

3-5576E 

01 

3 

• 2977E 

00 

5 

6* 

346 

5511 .345 

6 

• 8 469E 

8.2156E-05 

1* 

1 9oOE 

05 

2* 

8000E-03 

1 * 

56o4E 

00 

4 . 4347E 

01 

1 

• 5 6 ^ 4 E 

00 

6 

9* 

435 

21231 .354 

3 

•0420E-C4 

1 .369?E-0^ 

6 • 

2300E 

04 

7 • 

2000E-03 

9t 

1 659E 

-01 

i • 0 1 3 1 E 

01 

9 

• 1 659E- 

0 1 

7 

lO- 

ICO 

8344.861 

7 

• 55 1 4E-c*5 

1-1329E-0'* 

2* 

5375E 

0* 

8. 

0000E-03 

2* 

9290E 

00 

2-9136E 

01 

2 

•929oE 

00 

8 

ll* 

4 10 

20773-187 

2 

• 48o7E’-o4 

7.4421E-05 

3 • 

8500E 

05 

7. 

60C0E-03 

5* 

8 1 4SE 

00 

6 • o?S6E 

01 

*\ 

•8i*8E 

00 

9 

6* 

425 

1 1632-095 

1 

•2576E-03 

2 - 2644E”C^ 

1 * 

34Q5E 

05 

3 • 

4000E-03 

7. 

38 1 OE 

00 

1 .7275E 

0? 

7 

•38jnE 

00 

10 

9 * 

965 

5564.723 . 

1 

• 1 757E-04 

2*35 1 4E-C4 

3* 

5000E 

06 

4 • 

0000E-03 

4* 

7669E 

00 

9.4834E 

01 

4 

•7669E 

00 

11 

1 1 * 

40? 

21929.725 

• 2 

•4296E-Q4 

9.7147 E -o5 

4- 

4625E 

03 

9. 

0000E-03 

7 1 

6950E 

00 

6.8039E 

01 

7 

•6950£ 

00 

12 

6* 

693 

6218.612 

8 

•3724E-04 

2 •0094E-0'* 

1* 

3 1 42E 

05 

3. 

7000E-03 

3. 

7876E 

00 

8 . 1 4 6 1 F 

01 

3 

•7876F 

00 

13 

10* 

227 

5675.929 

1 

•2113E-04 

1 • 938oE "0^ 

1* 

0920E 

0* 

5. 

'♦C00E-03 

2 * 

2447E 

00 

3. 3o79E 

01 

2 

• 2 <♦ 'i 7 E 

00 

14 

11* 

535 

5920*581 

1 

•3441E-C4 

'♦•0330E-05 

2 • 

1000E 

0^ 

6 • 

2000E-03 

7. 

7 630E 1 

■oi 

9 • 9638E 

00 

7 

• 7 67qE- 

01 

15 

7. 

265 

124q6.o86 

1 

• 7 4 58E- q3 

7 . 3323E"0^ 

1* 

4385E 

0^ 

4* 

7000E-03 

2* 

6977 E 

03 

4 . 5 6 7 6 E 

02 

2 

• 6 9 7 7 E 

01 

1 6 

10* 

793 

S37l .55 q 

1 

• 8 1 24E »o4 

1 •8124E-0'* 

5. 

4o75e 

0 5 

1 • 

0000E-02 

5 f 

1 324e 

00 

^• 08 A 2 E 

01 

5 

• 1 3?4E 

00 

17 

11* 

776 

5164.383 

1 

•7734E-Q4 

1 •7734E-0'* 

5. 

8975E 

0^ 

4. 

9000E-03 

3. 

1260E 

00 

5 . 0 767E 

01 

3 

• 126oE 

00 

is 

4* 

858 

3251.649 

4 

•1564E-Q5 

2-1051E-03 

3. 

8150E 

OS 

5. 

7000E-03 

2. 

1 592E 

01 

3.Q703E 

02 

p 

•1992E 

01 

19 

8. 

6?S 

544 0 . 173 

1 

•0S43E-Q3 

7.2642E-0'* 

1 • 

0675E 

0 5 

6 • 

2000E-03 

1* 

4369E 

01 

1.8443E 

02 

1 

•4369E 

01 

20 

9- 

So 2 

17583.814 

6 

•2731E-04 

1 »88l5E-0^ 

1 • 

7675e 

05 

6 • 

7 OOOE -03 

1 • 

q437e 

01 

1.2396E 

02 

1 

*0'»37E 

01 

21 

11* 

576 

16987.752 

3 

• 3924E -q4 

2.2o52e-0^ 

2. 

3800E 

0^ 

1 ♦ 

5000E-02 

1 • 

3657e 

01 

7 . 2 4 5 1 E 

01 

1 

•3657E 

01 


vj 
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FORCE 
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XN 

CM3 


TEST NO# 7 SATURN V DT V CONFIG# II LONGITUDINAL F8RCE LEVEL (1967) 

ZETA ’ 


XE 

CM] 


genu mass 

CKG-SEC2/M] 


D C«PI*Xe*F] 
CN-M] 


T D I 

CN-M3 Cn-M] 


1 

6 

•220 

5675.929 

6 

•1290E-0" 

6.3240E-05 

7 

• 6825E 

04 

* OOOOE 

00 

1 

•1277E 

00 

2 

• 2039E 

01 

•OOOOS 

00 

2 

6 

• 22b 

4150.189 

4 

•6884E-04 

4 .7058E-05 

7 

• 6325E 

04 

• OOOOE 

00 

6 

• 1354E 

-01 

1 

•2917E 

01 

•OOOOE 

00 

3 

6 

• 230 

2302.379 

3 

• 1793E-0 4 

3 • 1756E-05 

7 

•6825e 

04 

•OOOOE 

00 

2 

•7958e 

-01 

5 

•9493E 

00 

•OOOOE 

00 

4 

9 

.180 

21636.142 

7 

.1529E-04 

1 .5542E-04 

9 

. 5900E 

04 

.OOOOE 

00 

1 

• 0564E 

01 

8 

. 162 IE 

01 

•OOOOE 

00 

5 

9 

.182 

16026.937 

5 

•0149E-04 

1 .1072E-04 

Q 

. 5900E 

04 

•OOOOE 

00 

5 

•5749E 

00 

4 

• 0137E 

01 

•OOOOE 

00 

6 

9 

.184 

10617.901 

3 

.723CE-04 

7.3955E-05 

9 

• 5900E 

04 

.OOOOE 

00 

2 

• 4 66 9E 

00 

2 

• 2 1 3 1 E 

01 

•OOOOE 

GO 

7 

11 

• 20(? 

21934.173 

3 

* 6969E-04 

9.891 7E-05 

1 

• 3177E 

05 

•OOOOE 

00 

6 

« 8 1 &2e 

00 

4 

• 4594E 

01 

•OOOOE 

00 

8 

11 

• 190 

16089.212 

3 

•0202E-04 

7.2544E-05 

1 

• 3 177E 

05 

•OOOOE 

00 

3 

•6668E 

00 

2 

• 9708E 

01 

•OOOOE 

00 

9 

1 1 

• 195 

10795.330 

2 

• 1 632E-04 

4.5786E-04 

1 

•3177E 

05 

•OOOOE 

00 

1 

•5529E 

01 

1 

•5255E 

01 

•OOOOE 

00 

10 

6 

• 343 

5876.099 

7 

• 2588E-0 4 

9.6403E-05 

1 

•1900E 

05 

•OOOOE 

00 

1 

•7796E 

00 

4 

•9796E 

01 

•oooo£ 

00 

11 

6 

• 3'.r, 

4052.323 

•5 

•0061E-0" 

6.5598E-05 

1 

• 1900E 

05 

•OOOOE 

00 

s 

•351 IE 

-01 

2 

• 3699E 

01 

•OOOOE 

DC 

12 

6 

.346 

2771.241 

3 

•4S24E-0 4 

4-5422E-05 

1 

•1900E 

05 

•OOOOE 

00 

3 

•9545E 

-01 

1 

• 1 472E 

01 

• 0 0 0 0 E 

00 

13 

9 

• 410 

21177.975 

2 

•7714E-0 4 

1 -4077E-04 

6 

♦2300E 

04 

•OOOOE 

00 

9 

* 3660E 

00 

8 

• 3636E 

00 

•OCTQE 

00 

14 

9 

-411 

15920*179 

1 

• 2979E-0 4 

9.4298E-05 

6 

•2300E 

04 

•OOOOE 

00 

4 

•7163E 

00 

1 

•834SE 

00 

•OOOOE 

00 

15 

9 

.412 

1 1076. 068 

7 

•0591E-05 

6 • 1 563E-05 

6 

•2300E 

04 

•OOOOE 

00 

2 

* 1422E 

00 

5 

• '•285E- 

-01 

•0000E 

00 

16 

1 0 

. C67 

8456*066 

7 

•3528E-05 

1 .3126E-04 

2 

. 5375E 

06 

•OOOOE 

00 

3 

• 4870E 

00 

2 

• 7 4 4 3 E 

01 

•OOOOE 

00 

17 

10 

• 068 

6276.438 

6 

•5457E-05 

1 .04S8E-04 

2 

• 5375E 

0^ 

•OOOOE 

00 

2 

• 0679E 

00 

2 

• 1 753E 

01 

•OOCOE 

00 

is 

10 

•080 

4225. 8 0 9 

4 

•5684E-05 

6.1733E-05 

2 

• 5375E 

0* 

•OOOOE 

00 

8 

• 1 956E 

-01 

1 

•0622E 

01 

*0000^ 

00 

19 

1 1 

.535 

21435.972 

2 

.4795E-0" 

7.5914E-05 

3 

• 85qoE 

0 5 

• OOOOE 

00 

5 

* 1 1 23E 

00 

6 

• 2166E 

01 

•ooonE 

00 

20 

1 1 

.437 

1 572q. 009 

1 

. 4805E-0 4 

5 r 61 73E-Q5 

3 

•8500E 

05 

•OOOOE 

00 

2 

•7742E 

00 

2 

• 1739E 

01 

•OOOOE 

00 

21 

6 

• 419 

10809*175 

1 

•1904E-03 

1 .9484E-04 

1 

. 34q5e 

05 

•OOOOE 

00 

6 

• 6 1 65e 

00 

1 

• 5451E 

0? 

•oonoE 

00 

22 

6 

.420 

8 042.382 

9 

•I444E-0 4 

1 .5059E-0* 

1 

• 34q5e 

05 

•OOOOE 

00 

3 

•80 4 8E 

00 

9 

• 1 1 96E 

01 

•OOOOE 

00 

23 

6 

.421 

5333.416 

6 

•2328E-0 4 

1 .Q258E-0 4 

1 

•3 4 o5E 

05 

•OOOOE 

00 

1 

* 7 1 8 8 E 

00 

4 

• 238 IE 

01 

•COOCE 

00 

24 

9 

.955 

6094.061 

1 

•2R26E-0 4 

2.8379E-C 4 

3 

•5000E 

0* 

•OOOOE 

00 

5 

■ 433 1 E 

00 

1 

• 1 263E 

02 

•cooo£ 

00 

25 

9 

•96q 

4377.048 

9 

•7968E-05 

2.0126E-0 4 

3 

•5000E 

0$ 

•OOOOE 

00 

2 

• 7675E 

00 

6 

• 5779E 

01 

•OOOCE 

00 

26 

9 

•961 

2966.963 

6 

•2253E-05 

1 .2759e-o4 

3 

•5000E 

0* 

•OOOOE 

00 

1 

•1893E 

00 

2 

•6565E 

05 

•OOOOE 

00 

27 

1 1 

• 53c 

15723.906 

1 

•3677E-Q4 

6.5437e-q5 

4 

• 4625E 

05 

•OOOOE 

00 

3 

• 2335E 

00 

2 

•1905E 

01 

. OOOOE 

00 

28 

11 

.525 

10684.624 

1 

.3467E-C" 

4.3753E-05 

4 

• 4625E 

05 

OOOOE 

00 

1 

.4687e 

00 

2 

• 1221E 

01 

•OOOOE 

00 

29 

6 

.7 0 5 

6458.815 

8 

•6407E-04 

2 .0089E-04 

1 

. 3142E 

05 

•OOOOE 

00 

4 

• 0763E 

00 

8 

•7Q76E 

01 

• 00 .OE 

00 

30 

6 

.707 

4910*835 

6 

•5780E-0" 

1 .5314 E -o4 

1 

• 3 1 4 2e 

05 

•OOOOE 

00 

2 

• 3626e 

00 

5 

* 0 4 96£ 

01 

•OOOOE 

00 

31 

6 

• 7 r 9 

3180* "77 

4 

.4984E-04 

1 *0571 E*0 4 

1 

«3142e 

05 

•OOOOE 

00 

1 

• 0562e 

00 

2 

• 3629E 

01 

•OOOOE 

00 

32 

1C 

*29S 

5925.029 

2 

•1334 E -q4 

2 . 896 1 E -O 4 

1 

• 0920E 

0^ 

• OOOOE 

00 

5 

•39q8e 

00 

1 

• 0 4 06E 

02 

■ ’COCOE 

00 

33 

to 

• 3lo 

2851. 3q9 

9 

•8926 E - 0 5 

1 .3125 e - 0 4 

1 
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TEST 

N0- 8 SATURN 

V OTV CONFIG, II CO MSFC 201 

PITCH + YAW FREQ. 

SWEEP (1967) 



fregJ^y 

FORCE 

XN 

XE 

gem mass 

^Pzeta 

D t = PI*XE 

*E] 

T 

DC* 

* I 

A*T} 

[HZ] 

CNJ 

CM3 

CM] 

CKG-SEC2/M] 


CN-M] 


CN-M] 


Cn:-v]* 


1 ♦ 6 &o 

4830-767 

2 - 8126 E-Q 2 

9.o94Se-q4 

5-4950E 

03 
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1 e 3371 E 

01 

2. 3645E 

oa 
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01 
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01 
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02 
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01 
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TEST 

N0 

8 SATURN 

V 

DTV CQMF 

JG. II Cfl 

MSFC 201 

PITCH + YAW FORCE 

LEVEL (1967) 

4*Pl*^j^ 



F0RCE 


XN 


XE 

gene mass 

^ZETA 

* 

DC 

=Pl*xr*FJ 


T 

D O 

\*T) 


CH23 

CM? 


CM] 


CM] 

CKG-SEC2/M] 




CN-M] 



CN-M] 


CN-M] 


i 
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5* 

4950E 

03 

•OOOOE 

00 

4 

•0139E 

00 

7* 

2122E 

01 

•OOOOE 
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1 

•5262E-04 

2* 

3625E 

03 

• OOOOE 

00 

7 

• 7q 1 6E 

00 

1 

1 525E 

02 

•oo'inE 

00 


5*908 

12005.746 

6 

•9514E-04 
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6 • 

1950E 

03 

•OOOOE 
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TEST 

NB, 9 SATURN. 

: v DTV C.6NF I G • ! I C9 MSFC 201 T0RSI0NAL FRE3* S^EFP (1967) 


4*Pl4^A»T3 

FRECL®^.Y 

T6RQUE 

XN 

XE 

GF.NL MASS 

^^ZETA ’ 0C = PI*XE*F] 

T 

oc* 

cHzr 

CN-M3 

CRaD] 

CRAD3 

CKG-M23 

CN-M] 

CN-M] 


[N-W 

7*893 

18099.396 

4.8000E-05 

4.3735E-06 

2 • 2A83E 05 
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1‘OOOOE-O 1 * 
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sn 
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test NO. 9 SATURN V OTV C6NFJG. II CO MSFC 201 TORSIONAL FORCE LEVEL { 1967 J 



freqi^P:y 

TORQUE 

XN 

XE 

genl mass ^Pzeta 

* 


DC*PI*XE*F3 

T DO 

4«Pl*^^A*n 


inrr 

CN-M] 

CRAO] 

[RADI 

CKG-M2I 



CN-M3 

CN-M] 

cn-mi 

1 

8 • 6q5 

1 7 A 46 . 394 

1 .0550E-0^ 

9 .5l7if-o6 

1.4122E c 5 • OOOOE 

00 

t 

5.2163E-01 

2. 2974 E 00 

•oocoe 00 

2 

8 • 63q 

12784.410 

9.1097E-05 

7.9960E-06 

1 • 4 1 22E 05 • OOOOE 

00 


3.2114E-01 
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•OOOOE 00 

3 

8.675 

8122.426 

6.1S14E-05 
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1 • 4 1 22E 05 • OOOOE 

00 
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8.015SE-01 

•OOOOE 00 
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2.6400E-0& 

4 .45o2p-o6 

7 . 0951 £ 07 .OOOOE 

00 
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1.3192E 00 

•OOOOE 00 

5 

1 1.632 

9323.968 

2.1187E-06 

3.1352E-C6 

7.0951E 07 .OOOOE 

00 


9.3301E-02 

8.5C58E-01 

•OOOOE 00 

6 

11.655 

4661.984 

1.4441E-Q6 
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7 • 0951 E 07 .OOOOE 

00 
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•OOOOE 00 
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TEST N6. 10 SATURN V DTV CONFIG. II C3 MSFC 201 LONGITUDINAL FREGU S^EEP (1967) 



fred^^:y 

FORCE 

XN 

XE 

genl mass 

^FZETA 


D C = PI*XE 

*n 

T 

D C c 

4*Pl«!^ 

A * T] 


CHZ3 

Cn3 

Cm] 

CM3 

CKG-SEC2/M3 



CN-M3 


CN-M3 


Cn-m] 


1 

6.211 

5667.032 

6.li75E-0^ 

7.3423E-05 

9.1175E 

0<* 

3.4000E-03 


1-1101E 

00 

2.5982E 

01 

1 * 1 10lE 

00 

a 

9-195 

19802.A75 

5#6703E-0'+ 

8.5021E-05 

1.2197E 

05 

8.9000E-03 


7.3199E 

00 

6.5450E 

01 

7*3l"9E 

00 

3 

9*758 

17183.474 

2 • 8702E -q5 

2*0055E-04 

4-6200E 

07 

l'lOOOE-02 


9.8879E 

00 

Z * 1 532E 

01 

9 * 8879E 

00 

4 

11.212 

17232. 4 0 4 

3 • 3&03E -Q4- 

8*1145e-o5 

1 . A455E 

05 

6.8000E-03 


3. 2580E 

00 

3-S123E 

01 

3*25«oF 

00 

5 

8*628 

4790*733 

9-2433E-0^ 

6 • lSl 4E-C-4 

1 • C097E 

05 

6-7000E-03 


1 * 0873E 

01 

1.2677E 

0? 

1 ‘0673E 

01 

6 

9.513 

17232.404 

6*01 13E-0 4 * 

1 .5632E-04 

2 • 2575E 


6.3000E-03 


1 . 1536E 

01 

1.4572E 

02 

1 * 1536E 

01 

7 

11.587 

17499.297 

3 • 2376E- 0^ 

2.431 1 E*Q4 

2-5200E 

05 

1.5000E-02 


1 .3195E 

01 

7 • 0002E 

01 

1 -3195E 

01 


6 ^ 



TEST N0 • 10 SATURN V DTV CONFIG. II Cp MSFC 201 LONGITUDINAL FORCE LEVEL (1967) 



FREG'!^Y 

FORCE 

XN 

XE 

genl mass 

• z £ x. 

* 


DC=PI*XE 

*F3 

T 

DC* 

4*P W^F 

A # T3 


CHZ1 

CN3 

CM] 

CM] 

CKG-SEC2/M] 



4 

CN-M] 


CN-M] 


Cn-m] 


1 

6*212 

5168.832 

5.6AS2E-0'* 

6.2279E-05 

9-1175E 

04 

•OOOOE 

00 


1.0U3E 

00 

2.2156E 

01 

•OOOOE 

00 

2 

6*211 

3932.226 

4.9296E-0A 

5.2523E-C5 

9* 1 175E 

04 

• OOOOE 

00 


6.4S84E 

-01 

1 • 687 1 E 

01 

•0000- 

00 

3 

9*170 

1 7445.919 

4.U1EE-0 4 

1 . 1549^-04 

1*2197 E 

05 

•OOOOE 

00 


6 * 3297;r 

00 

3.4230E 

01 

•'OOOOE 

00 

4 

9*169 

12837.563 

2*87S2r-Q4 

8*1257 E -o5 

l*2l97f 

05 

•OOOOE 

00 


3 • 277 l [t 

00 

1 • 6 7 7 6 E 

01 

•OOOOE 

00 

5 

9-168 

8718.511 

1 .8693E-0A 

5 .4421 E-05 

1 *2197e 

05 

.OOOOE 

00 


1 . 4906E 

00 

7.071 IE 

00 

•OOOOE 

00 

6 

1 1.200 

17343.610 

3.1718E-04 

6 .9S99E-05 

1 *4455 E 

05 

.OOOOE 

00 


3,80S6e 

00 

3.6C08E 

01 

♦OOOOE 

00 

7 

11*199 

12579.566 

2.4397E-0'* 

5.23l5 E -05 

1 . 4455 E 

05 

.OOOOE 

00 


2.0^>75e 

00 

2.21S3E 

01 

•OOOOE 

00 

S 

11.190 

8340.412 

1 .6448E-0A 

3 .5509E-CT5 

1 .44P.5E 

05 

•OOOOE 

00 


9.3338E' 

*01 

3 • 6 66 1 E 

00 

♦OOOOE 

00 

9 

8*628 

^ 6 7 0 • 6 3 1 

9*3?70E-04 

6*1 070E-04 

1.0097E 

05 

•OOOOE 

00 


8.9609E 

DO 

1.3102E 

02 

•00 n 0E 

00 

10 

8*626 

4261.395 

8 » 5584E-04 

5.5777E-C4 

1.0097E 

05 

.OOOOE 

00 


7.4672E 

00 

1.0863E 

02 

•OOOOE 

00 

1 1 

8 * 6C 4 

2579.968 

' 4.8468E-04 

3*1731 E"04 

1 • 0097E 

05 

•OOOOE 

00 


2 * 571 9e 

00 

3.4S24E 

01 

♦OOOOE 

00 

12 

9.553 

17116.751 

5.0726E-0* 

2.6793E-0'* 

2-2575E 

05 

•OOOOE 

00 


1 *44o7E 

01 

1 . 0 A 64 E 

02 

• cooqE 

00 

13 

9-549 

13015.492 

3-9792E-0'* 

2*l793E-0“ 

2*2575 E 

0 5 

•OOOOE 

00 


8.91 10E 

00 

6.4337E 

01 

• CO 

00 

1* 

9-544 

8553.927 

2.6777E-QA 

1 • 45l 1 E~0* 

2*2575E 

05 

•OOOOE 

00 


3.S995E 

00 

2.9103E 

01 

* 0 0 0 0 E 

00 

l5 

1 1 -634 

1.7281 *335 

3.0873E-0 4 

1 *9343e-0^ 

2*52oOE 

05 

•OOOOE 

00 


1 *0773 E 

01 

6 . 4 1 73E 

03 

•OOOOE 

00 

16 

11 .645 

13300* 1 78 

2.8487E-0'* 

1 *5339 E -o4 

2 • 52Q0E 

05 

•OOOOE 

00 


6.4926 E 

00 

5. 4738E 

01 

•oooqe 

00 

17 

1 1*635 

8727. 4q8 

1.8844E- 0 4 

9.7314E-05 

2 • 52Q0E 

05 

•OOOOE 

00 


2 • 668 1 E 

00 

2 • 39 1 oE 

01 

•OOOOE 

00 




TEST N0 • 11 S-IV-B-O/ LEM, AND APfiLLO PITCH * YAW (1966) 



FREOl^Jtv 

FORCE 

XN 

XE 

genl mass 

^ZETA 


DC=PI*XE*F} 

T 

D C 

4#PNWa#T] 


CHZJ^ 

CNl 

Cm} 

CM} 

CKQ-SEC2/M} 




CN-M} 

CN-M} 


Cm- v ] 

1 

4*980 

3229.000 

1 » 30 1 OE-03 

2.2546E-0'* 

1 • o555E 

04 

8 

•9000E-03 

i 

9.7814E-01 

8 . 7458E 

00 

9*7 81 4 E-01 

2 

9*660 

348i .000 

3 • 8500E-0 1 * 

9.0128E-05 

1 .6048E 

04 

1 

•1S00E-02 


6.4971E-01 

^ • 38 1 5E 

00 

6*49715-01 

3 

10*000 

3346.000 

2.9000E-0 1 * 

9.0973E-05 

1-9421E 

04 

1 

•2300E-02 


4.9832E-01 

3.224QE 

00 

4 *.98025-01 

4 

10*580 

3449.000 

2 » 2800E-04 

9.1314E-05 

3.7114E 

04 

1 

• 1400E-02 


6.1069E-01 

4 « 2 6 2 9 E 

00 

6*10695-01 

5 

14*800 

3305.000 

1.9000E-04 

1 •2969E-0 i » 

5.6828E 

04 

1 

•3000E-02 


1-4490E 00 

8.8700E 

00 

1 ♦ 4 4 9qE 00 

6 

5*200 

3705.000 

2-2820E-03 

1 .5561E-03 

1* 1824E 

04 

1 

•3700E-02 


5 . 6580E 00 

3 • 2S65E 

01 

5 * 65 9 oF CO 

7 

9*680 

3774.000 

3.9300E-0'* 

8.6578E-05 

1 • 3646E 

04 

9 

• 4000E-03 


4.60't8E-01 

3 » 8 9 8 3 E 

00 

4 • 6 04 85 - oi 

8 

1 0 * 9 1 0 

3736.000 

1 •6300E-0'* 

7 • 2568 E-o5 

2 • 2986E 

04 

1 

• 1500E-02 


2-0736E-01 

1 -4349E 

00 

2*C716F-01 

9 

15*000 

3373.000 

2.0000E-0'* 

1 .1468E-0'* 

2 • 2630E 

04 

1 

•OSOOE-02 


5.4562E-01 

* * 0~ J3E 

00 

5 *45A2E-0l 

10 

4*900 

3371.000 . 

1 .9850E-03 

3 •5015E-C'* 

1-0284E 

04 

1 

•9400E-02 


4.6S18E 00 

1 -9205E 

01 

4 • 6 8 } 8 E 00 

11 

9*350 

3609.000 

• 3*61 OOE-O^ 

9.9347E-05 

2.3731E 

Q4 

1 

•0900E-02 


7.3727E-01 

5.3826E 

00 

7*37275-01 

12 

1 O’-CO 

3395.000 

2.3100E-0'* 

6 .2693E-05 

3.Q539E 

04 

7 

•6C00E-03 


3.6499E-01 

3 . 8 2 1 7 E 

00 

3*64995-01 

13 

14*420 

3502-000 

1.590CE-0'* 

1 . 1 64 4 e -q4 

1 .8682E 

05 

9 

•5000E-03 


2 * 3 1 ^ 2E 00 

1.9385E 

01 

2 • 3 1 fl 2 5 00 


5 • 09q 

3346.000 

8.310CE-04 

7 •560 / *E"0 < * 

1-2523E 

04 

1 

•4400E-02 


8.Q029E-01 

* * 4226E 

00 

8*00*95-01 

15 

9*440 

3415.000 

5.4600E-0'* 

1 •3426E-0'* 

1 -0630E 

04 

9 

• OOOOE -03 


6.3C-44E-01 

5.5743E 

00 

6 «3044E-01 

16 

10*950 

3546.000 

3.5700E-0'*- 

9.2570E-05 

1 .4796E 

04 

1 

• 1400E-02 


6.3937E-01 

^.4631E 

00 

6 • 3937E-0 1 

17 

13*490 

3732.000 

1*51 OOE-O't 

6.909Se^05 

5.5456E 

05 

1 

• OOOOE - 02 


5.707SE 00 

4.542 IE 

01 

5*70785 00 

13 

1 5 * 2 0 0 

3907.000 

1 ♦5800E-0'* 

4 .7447E-C 5 

2-1740E 

04 

1 

• 0°00E - 02 


3.39Q2E-01 

2 • 4 75 1 E 

00 

3*39-25-01 


S' 



^ in cr> 


freqJ^y 


TEST 

N0. 12 S-IV 

-0-0/ LEM/ 

AND AP6LL8 T0RS I 0NAI. 

( 1966) 




TORQUE 

XN 

XE 

GENL MASS 

^ZETA 

DC«PI*XE*FJ 

T 

D[=4*PI#^Pa»T] 

CHZ1 

CN-M] 

CRaDJ 

CRADl 

CKG-M2] 


cn-mJ 

CN-M1 


CN-mJ 

4 • 90Q 

16341.000 

8-6926E-04 

3*2250E-05 

2*9000E o'* 

l*3600E-02 

1 « 6556E 00 

1 * 0385E 

01 

1 • 7 7 a 9 E 00 

5*610 

16330*000 

l*640lE-o3 

1 -6237E-05 

1 • 5340E O'* 

3*0600E-02 

8.3297E-01 

2 • 5633E 

01 

9*8566E 00 

9.330 

15750.000 

9.3619E-05 

5 . 1743E-05 

3.CS61E 06 

1 .4600E-02 

2.5603E 00 

4 * 6477£ 

01 

s *527oE 00 

5*160 

17027.000 

9.039SE-04 

1 . 5625E -04 

1 • 07C4E 05 

2-5900E-02 

8 * 3 5 8 3 E 00 

4 • 643 1 E 

01 

1*51=8E 01 

5*700 

16999.000 

1 .6913E-03 

9.8532E-05 

9-9050E 03 

1 .7900E-02 

5.2620E 00 

1 *6905E 

01 

3 • 8C-R6E 00 

9*430 

16279-000 

7.<t'*03E-05 

3.8348E-05 

3 • 8921 E 06 

•OOOOE 00 

1 » 96 1 2E 00 

3.7S20E 

01 

*0000- 00 


ss 

N 



- 

• 

— 

TEST 

N9. 13 SAD- 

202 UPPER 

ST 

PITCH ♦ YAW (1966) 




m 


FRECUENCY 

FORCE 

XN 

XE 

GENL f*ASS 

ZETA 

OC*?I*XE->n 

T 

0 

Ta.T) 


tHZJ 

Cn3 

I Ml 

tM] 

CKG-SEC2/ M 3 


. CN-H3 1 

CN-R) 



1 

2*090 

3022*000 

1.087CE-03 

4.4784E-0'* 

9.7334E 

05 

3* 1C00E-02 

3.8629E 01 

9*91622 

01 

3*86295 ol 

2 

7*810 

5199.000 

2*0590E-03 

7.3280E-04 

5.0200E 

04 

1 •0 < *00E-02 

3.3“88E 01 

2*562^2 

02 

3*34»3E 01 

3 

1 0 * 090 

5466*000 

4.5000E-06 

1 .81<*2E-05 

7.5717E 

06 

8.6000E-03 

3* 32992-02 

3 • 08 1 32 ' 

“01 

3 • 

4 

1 C* 820 

5482*000 

1.4100E-04 

6.7116E-05 

1 *5866E 

05 

1 *68005-02 

1*53892 00 

7 • 25932 

00 

l*53*9E oo 

5 

14*310 

4478*000 

1 .4000E-0 1 * 

8.67302-05 

2«8320E 

04 

1.4CC0E-O2 

4.Q1692-01 

2.25332 

oc 

4*01692-01 

6 

8*160 

2327*000 

1 «645oE-o3 

4.5912E-0'* 

2 * 11642 

04 

9*8C00E-03 

9.26992, 00 

7.52732 

01 

9*269=5 CO 

7 

14*0^0 ' 

2380*000 

l*520CE-0‘‘ 

5.3384E-05 

1*9634E 

04 

S*5CC0E*03 

1*393‘*2-01 

1*77262 

CO 

1*89342-01 

8 

8*49o 

1988*000 

1.593CE-03 

7.5359E-04 

1*35i7E 

04 

l*03C0E-02 

6 * 3 l 6 S£ 33 

4 * 3So4E 

01 

6*3ia 92 oo 

9 

14*270 

2277.000 

1.070CE-0'* 

3.8456E-C5 

1.5436E 

04 

1.4QC0E-02 

1.24972-01 

7*10362 

*01 

1*24972-01 

1 C 

2*030 

3345*000 

8 * 6400E-C4 

3*1 1 13E-C4 

1.1934E 

06 

1 . 8000 E -02 

1.72092 01 

7*60822 

01 

l*72o e 5 01 

11 

7*690 

3536.000 

2.0340E-03 

5 ,2233E”0 < * 

4.7160E 

04 

8.40C0E-03 

2.4C41E 01 

2*27752 

02 

2*4o 1 *lE Cl 

-12 

10*380 

3782*000 

2.4Q00E-05 

1 *2057E-0^ 

2*1985E 

07 

7.2D00E-03 

2.4367E 00 

2.69322 

01 

2*43672 00 

13 

11*010 

3533.000 

1 . 2600 E- 0 *' 

4.6696E-05 

4.57q9E 

05 

1 *02C0E-02 

2.22562 00 

1 .726i*2 

01 

2*22562 CO 

14 

8*060 

2303.000 

1.8490E-03 

4,i787E-0^ 

2*17c5E 

04 

1 .0?00E-02 

1*29142 01 

9«5l55E 

01 

1 • 2 9 1 * E 01 

1 5 

13*580 

2463.000 

4.4Q00E-05 

3.3546E-05 

1 . 1366E 

06 

l.OSCOE -02 

1*q871E 00 

8 . o 1022 

00 

1 *0871 E 00 

16 

8*310 

2014.000 

2.1660E-Q3 

8.2005E-04 

1-2720E 

04 

9.9000E-03 

1.0120E 01 

8.13462 

01 

1*C1205 Cl 


00 


‘.DOO'-joai^ujro*-* 



CHZ3 

7*850 
4-990 
10*800 
5*110 
8*190 
10*800 
5 * 1*0 
8 • 49 q 
1 1 *200 


TORQUE 

tN-MJ 

10564 *000 
10 * 36*000 

10864.000 
11430*000 

11268.000 

10963.000 

12980.000 

12042.000 

12149.000 



S ' 1 

'■o 






TEST 

N0 ♦ 15 3eE I N 

•'3 747 ( 1969 


FReJ^N'CY 

FORCE 

XN 

XE 

GENU MASS ^ 


[HZ] 

CN3 

CM] 

CM] 

CKG-SEC2/M] 

1 

1*015 

2499*000 

7 • 0500E-05 

2.4964E-05 

4 • 5587E 

Q5 

2 

1*635 

2499*000 

7.0500E-05 

2. 4964^-05 

4.5587E 

05 

3 

2*100 

2499.000 

7 • 0500E -05 

2 .4964E-05 

4*55g7E 

Oo 

4 

2*^25 

2499*000 

7.0500E-05 

2.4964E-05 

4.5587E 

05 

5 

2*500 

2499*000 

7«0500E-c5 

2.4964E-05 

4.5587E 

03 

6 

3*16o 

2499*000 

7 - 0500E-05 

2.4964E-05 

4 ♦ 5 5 8 7 E 

05 

7 

3*46o 

2499*000 

7 • 0500E - 05 

2.4964E705 

4.5587E 

05 

8 

*948 

2499.000 

7 »0500E-05 

2 *4964E"05 

4 • 5 5 8 7 E 

05 

9 

1.775 

2499.000 

7.0500E-05 

2 ,4964£-o5 

4 * 5587r 

05 

10 

1-900 

2499*000 

7 • OSOCE -05 

2*4964 E -05 

4.5587E 

05 

11 

2*<i90 

2499*000 

7 • 0500E - 05 

2 • 4964 E -05 

4.5587E 

05 

12 

3*090 

2499*000 

7 • 05C0E - q5 

2.4964E-05 

4.5587E 

05 

13 

3*6?o 

2499*000 

7 • 0500E-05 

2-4964E-05 

4 • 5 5 8 7 E 

q5 

14 

4 • 820 

2499*000 

7 *05002-05 

2 .4964E-05 

4 • 558 7 E 

05 

15 

5 • 060 

2499*000 

7 • 050CE - 05 

2 *4964f -o5 

4.5587e 

05 

16 

5*460 

2499*000 

7 • 0500E -05 

2*4964E-*05 

4 • 5587E 

05 


o 





APPENDIX B 


DIGITAL COMPUTER PROGRAMS 



C6MM0N D , T , W , X 

DIMENSION 0(3‘tC),T(3'!0),W(340),X(3W) 
DIMENSION Z<4),ARRAY(14) 

NMAX = 340 

READ 5, INDEX, 12, IS 


INDEX 

INDEX 

INDEX 

INDEX 


1 - D=C*X1**A1 

2 - D=C*<X1**A1 

4 - D=A1+A2»X1 ■ 

5 - 0 = C« (M**A; ) 


* (X2*»A2) 

A3<X1»X1 
► ( W * * A 2 ) * 


(X**A3> (UNLESS 12-7) 


:*t**ai 
C*X# *A 1 
C* W»*A 1 
C»M**A1 

Cl ( M * * A 1 ) * ( X**A2 ) 

C» <M*« A1 ) *( W**A2) 
C»(X*»A1)*(W»*A2) 
C*(T**A1M(X/L)**A2 SR C * ( M** A 1 ) 

5, (APRAY(I), I = 1,14),IST9P 


* (W»*A2) * ( X/l ) **A3 


1 READ 15, (APRAY(I), I = 1 , 1 4 > , I ST9P 
15 FORMAT ( 13AA,A1, ID 
IF < ISTSP-9)?5, 101,25 

25 PRINT 110, ( ARRAY ( I ) , I = 1,14), INDEX, 12 

110 FORMAT ( 1 M 1 , 35x, 13A6,A2/ 55X, 7HINDEX = , 12, 5H 12 =, I?/ 

* 55X, 1 OH I NPUT DATA// 6X, 9HFRE3UENCY, 5X,5HFSRCE, 9X, 

•» 2HXN, 10X , 2'" J XE , 6X/9HGENL MASS, 7X, 2HX0» 8X, 4HZETA, 12X, 1HD, 1 1 X, 

« 1HT, 10X, 2HD2) 

I = 0 

4 READ 3, 1 1, IDUM, 13 

5 FSRMA T (1015) 


- INPUT UNITS - CPS, KSF, 


M, XE/XO, KG-S2/M 
l/LB, G/LB, LB-S2/IN 
M, XE/XO, KG 


LB, O/LB, G/LB, 
N, M, XE/XO, 
D, T IN N-M 
N-M, DEG, XE/XO, 


11=5 - CPS, N-M, DEG, XE/XO, KG-M2 

11=6 - CPS, IN-KIP, RAD, RAD/IN-LB, LB-S2-IN 

11=7 - CPS, LB, G/LB, G/LB, L3-S2-IN, R(IN),RE(IN) 

13=1 - D = P I •* XE* E 
13 = 2 - 0= 4*P I »ZET A*T 
13=4 - , GM = 1, D=P I »XE»F 

126 IF CI1-1) 132,112,112 
112 READ 10/CLGTH 

IF (19) 102,102,106 

102 IF (I1-?) 103,148,104 

103 PRINT 115 

115 FSRMAT ( 3x, 5 h < CPS ) , 7X, 5H ( KGF ) , 9X, 3h(M>, 6X, 7H(XE/X0i 

► 11M(M0-SEC2/M),30X,5H(N-M), 7X, 5H ( M ) , 8X,5H( N-M ) / ) 

GST9 106 

104 IF ( 11 - 4 ) 109,109,146 





35 


= 

3 6 

r 

11 = 

1 

• 

IKP'JT U- 

ITS - CPS> 

KG r / 

Mi XF/XO/ KG 

-S2/M 

= 

— ~t 
*. / 

r 

11 = 

2 

- 


C^Si 

LB/ 

G/L3i 3/LSi L5- 

S2/IN 

= 

26 

C 

I 3. = 

3 

- 


c*s, 

N'i 

Mi XF/XQi 

KG 


39 

c 

1 

4 

- 


CF$, 

DiT 

IN N-M 



c 


11 = 

2 

- 


CP 3 / 

N->1/ 

DEG/ XF/XOi 

KG-^2 

= 

a: 

c 

i 1 = 

6 

• 


C~Si 

I r. ”< 1 P, "<aD, RAD/ 1 N-L 

3/ LB-S2-IN 

= 

u2 


11 = 

7 

- 


CP Si 

LB/ 

G/LBi G/LBi LB-S2 

-IN', R C I Ml , RE Cl 51 

i 

4 3 

c 









= 

u k 

r 

12 = 

i 

- 

0 = PJ *XE* 

F 




= 

it 0 

C 

12 = 

2 

- 

•jsU*PJ *7 

E T A x T 




= 

•* 6 

r 

T "< = 

4 

- 


P = P I *xr* F 




= 

k 7 

r 








i 

5 

US 

- 



I F 

CICU^-93 1 

?6i 1 01 > 126 



! 

= 

uV 


126 


IF 

[13-13 l 

3?, 112 , 11 ? 





PROGRAM TO FIND BEST -FIT EMPIRICAL DAMPING LAW 
THROUGH LEAST -SQUARES METHOD 




145 FORMAT ( 8X,5H(CPS), 8X, 3H(N), 10X, 3H(M), 6X, 7H(XE/X0), 7X, 

* 4H(KG> ,33X,5H(N^M.) ,7X,5H(N-M) , ?X , 5H ( N-M ) / ) 

GOTO 106 

146 IF ( 11-6)144,149/148 
144 PR I NT 147 

147 FORMAT ( 3X,5H(CPS), 7X, 5H1N-M), 8X, 5H t DEG ) , 5X, 7H(XE/X0), 6X, 

* 7H(KG-N2),31X, 5H<N-M), jt, 5H(N-M), 8X,5H(N-M)/> 

GOTO 106 

149 PRINT 151 • 

151 FORMAT ( 3X, 5H(CPS)» 6X, 8HUN-KIP) 6X,19H(RAD) (RAD/IN-lB), 2X 

* 12H(LB-S£C2-IN),23X, 5H(N-M), 7X, 5H(N-M), 8X, 5H(N-M)/) 

GOTO '"06 

148 PRINT 125 

125 FORMAT 8X / 5H ( CPS ) , 7X, aH ( L3 ) , 8X/ 6H t G/L3 ) , 6X, 6H < G/L3 ) , 3X, 

* 12H1LB-SECP/IN) , 29X,5H(N-M), 7X, SH ( N-M ) , 8*/ 5H(N-M)/) 

106 JNASS = IMaSS 

106 READ 10/ FRE3/XN/GM/F/XE/X0/ZETA 
10 FORMAT ( 7E10 • 3 ) 

IF (ABS(FREO) - 1 • E- 1 0 ) 4,4/128 
128 1*1+1 

IF t I -N w AX ) 129,129,130 
129 w( 1 ) * 6*2831 8»FREQ 
W2 = w ( I ) *W ( I ) 

TS = 9.SC665/W2 
IF (11-2) 107,57,54 
> 107 GM 1 = 9 • 80665*GM 

FI = 9 •S0665*F 
^ XNl = XN 

xei = xni*xe 

GOTO 43' 

54 I F t 1 1*- >58,59,55 

55 IF(H-6) 52,51,48 

52 XNl = XN«. 017453293 
GOTO 53- 

58 XNl = XN 

53 GM 1 = GM 
FI = F 

XEI = XN1*XE 
GOTO 43 

51 GM 1 = *11 298*GM 
FI = 1 12 *98 »F 
XNl * XN 

XEI = XE* F# 1 000 • 

GOTO 43 

48 GM1 = . 1129?.*GM 
FI = 24 . C3q? 46*F 
R = XC* • 0254 
R1=ZE.TA« .0254 
ZET A = 0 » 

XNl = TS*F*XN/R 
XEI * T 3 *F*XE/R 1 
GOTO 43 

57 GM1 = 1 75 • »GM 



FI = 4.44322«F 
XN1 = TS*F*XN 
XE1 = TS*F*XE 

43 D( I ) = 3 • 1 41 59*F 1 « XE 1 
TBAR= •5«'.s'2*GM1*XnI*XN1 
IF ( 12-1)46,44,70 

44 w2 = 1 • 

FRE0*1. 

G M 1 = 1 • 

GM =1 . 

GO TS 46 

70 IF (12 -3)71/72*73 

71 XNl = i • 

X N! = 1 -• 

GM 1 = 1* 

GM =1. 

G0 TS 46 

72 W2 «1. 

FR£3 =1. 

X \ 1 = 1 . 

XN = .!• 

G0 TS 46 

73 IF ( 12-5)74,75/76 

74 W ( I ) = XN1 
X N 1 = 1 • 

FREQ =1. 

G0 TS 46 

75 XN 1 = 1 • 

XN = 1 • 

GO T 0 46 

76 IF ( 12-7)77/46/46 

77 GM1=1. 

GM *1 • 

XNl=bQR T ( X NJ 1 ) 

46 XM • =• GM1*XN1*XN1 

T { I ) = ♦5»w2»XM 
GGTS 56 
59 T ( I ) = GM 
0(1) = F 

XM = 2«*GM/W2 

GM = 0* 

F • 0* 

56 D2 = 12*56636*2ETA*TBAR 
D I F F = ( D ( 1 ) -02 ) /P2» 100 • 

IF' ( 13-2) 117,118/117 

118 D ( 1 ) = 02 

IF ( ABS ( D ( I ) ) - l.E-10) 119,119/117 

119 I = I - 1 
GST0 106 

117 IF (IS) 116.116/1/4 

116 PRIM 120, I ,FREQ/F/XM/Xe,GM,XO/ZETA,0( I )/ T< 1 ) , D2 

120 FORMAT ( I4,F9.3,F13,3/ 1PE14,4,4e12.4/E14,4,E12.4/ 1X/E12.4) 
114 IF ( INDEX -4 ) 121,106/ 121 

121 D ( I ) = aL0G(D( I)) 




T ( I ) = ALSGtTt I ) ) 

W( I ) = ALBG (W( I ) ) 

X( I )=AL0G(XN1> 

IF ( INDEX - 2) 106* 60/122 
122 IF < INDEX-*) 106/ 106*170 

60 IF (12-7)61,65*65 

61 T < I >*ALB3<X>M 
G8 T8 106 

170 T ( I >=AL9G(GM1) 

I F ( 12-7)106,137,137 
137 X ( I ) = ALOG ( XNl/CLG~H) • 

GO T9 106 

65 W( I ) - '-LOG ( XN1/CLGTH) 

G9 TO 106 
132 NPTS = I 

135 CALL LESTSQiZ, NPTS, INDEX) 

GOTO 1 

130 PRINT 165, I 

165 FORMAT (//55X, 21HTO0 MANY DATA POINTS./ U> 
101 STOP 
END 





SUBROUTINE LESTSO (Z, N, InD) 

LINEAR LEAST SQUARES FIT OVeR TWO INDEPENDENT VAR I A3LES 
COMMON D,Xl,X2,X3 

DIMENSION D ( 200 ) , Xl(200), X2(200), X3 ( 200 ) , A ( A, 4 ) / 3(4,4), C(4,4) 
DIMENSION Z(4), S ( 4 ) » Y(4), <T<4), TOL ( 4 ) , AKT ( 4 ) 

CALCULATE COeFF ICI ENTS OF EQUATIONS 

DOS I = 1 , 4 
Y( I ) - 3. 

Z( I) = 0* 

s ( i ; = o . 

DO 5 J=!M 
C( I, J) = 0. 

5 A( I, J) = 0. 

A ( 1 , 1 ) = N 
DO 10 I = 1 , N 
X1SG = X 1 ( I ) * X 1 ( I ) 

A ( 1,2) = A( 1 ,2) + XI (I) 

A( 2,2) = A ( 2, 2 ) + X1SQ 
Y ( 1 ) = V { 1 ) + D ( I ) 

Y( 2) = Y ( 2 ) + D ( I > *X1 ( I > 

IF (IND-3) 7,7,8 

7 A < 1,3) = A( 1,3) + X2( I ) 

A ( 2 » 3 ) = A ( ?, 3 ) + Xl( I ) *X2 ( I ) 

A ( 3, 3 ) = A ( 3, 3 ) + X2 ( I ) «X2 ( I ) 

Y ( 3 ) ~ Y ( 3 ) + 0 ( I '*X2( I ) 

IF C IVD-E) 10,40,40 
40 At 1, 4 ) A( 1, 4 ) +X3( I ) 

A ( 2 , 4 ) * A ( 2, 4 ) +X1 ( I )*X3( I ) 

A(3,4)*a(3,4)+X2( I )*X3( I ) 

A ( 4>4 ) « A ( 4 , 4 ) + X 3 ( I )*X3( I ) 

Y ( 4 ) = Y ( 4 ) +D ( I ) *X3( I) 

GOTO 10 

8 I F ( IND-5)9,7,7 

9 X2( I lsXISQ 

A { 2,3) = A ( 2, 3 ) + X 1SG*X 1 ( I ) 

A ( 3, 3 ) = A ( 3, 3 ) + X 1SQ»X 1 SQ 
Y ( 3 ) = Y ( 3 ) + D( I ) *X 1 S3 

10 CONTINUE 

A (2,1) = A (1,2) 

At 3,2) = A( ?, 3) 

A ( 4 , 1 ) = A ( 1 , 4 ) 

A ( 4, 2 ) = A ( 2, 4 ) 

A { 4 , 3 ) = A ( 3 , 4 ) 

LIM = IND + 1 
IF (IND-3) 13,11. i2 

11 LIM = 2 
GOTO )3 

12 IF ( IND -3)14, 16, 16 
14 A( l»3)®/>(2,2) 



non 



GO TO 13 
. 16 LlM = 4 

13 A<3,1> = A(l, 3) 

PRINT 100, ( (A( I, J), J* 1/1), Y ( J ) / I = l/LIM) 

100 FORMAT ( ///A6X/3SHC6EFFTCIENTS OF LEAST SQUARES EQUATION// 

* 1PE36 • A/ E7u.4/E36.4,E16.4,E58.A/E36.4,2E16.4,E42.4/ 

* E36.A,3El6.A,E26.A) 

INVERT EQUATIONS . 

DO 140 2 = 5/4 
00 l-/' J = 1/4 
DO l'-C J = 1/3 
140 B( I » J) « A ( T , J ) 

IF (IND-2) 160/170/165 

165 IF( IND - 4 ) 160, 170/191 
170 CALL MATINV(A/3/3/0/DET,3) 

GOTO 18C 

160 CALL M A T I N V ( A ,2 / B / 0 / DE T / 2 ) 

GO IS ISO 

191 CALL MaJINVIA/ 4, B, 0, DET., 4) 

1«0 CONTINUE 

D6 20 I * l/LIM 
DO 20 K = l/LIM 
20 2 C I ) = Z.'I) + A(I/K)*Y(K) 

1 E « EXP(Z( 1 ) ) 

PRINT 120/ 7(1), E/ (Z<I), J = 2/LIM) 

120 FORMAT ( ///47X, 36HS0LUTI0N FOR LEAST SQUAP-ES CONSTANTS// 

« 55X/3HC = , 1PE16-4, 10X, 8HEXP(C) ■/ E16.4/ 

* . . ' j 4 X / 4 H A 1 =/E16.4/54v,LHA2 = / E 1 6 • 4/54 X, 4HA3 «/E16«4) 

PRINT "30 

130 FORMAT (///50X/30HERR0R OF LEAST SQUARES APPROX.) 

C 130 FORMAT {///50X.30HFRR0R OF LEAST SQUARES APPROX . //53X, 8H TERM NO., 
C « 8X, dHEPROR, 5X/ 16HL0G D3 AR/DaPPROX, 5X, 15HL0G D3AR/DCHANG/ ) 

DO 190 1=1/4 
HI * I 
KT( I )=0 

190 T 9l ( I )=HI«. 23026 
DO 25 I = l/N 

DEL = D ( I ) - Z(l) - Z(2)*X1(I) - Z ( 3 ) *X2 ( I ) - Z(4)#X3(I) 

ERR = EXP (-DEL)- i. 

DO 260 L = 1 / 4 
J= -L+5 

IF’ (ABS(DEL) -TCL(vJ) >210/210*260 
210 K T ( J ) = K I (J) + l 
260 CONTINUE 
25 CONTINUE 

. C 25 PRINT 131/ I.-ERR/ OEL, DEL2 
J 3 1 FORMAT ( I58/2F16. ■./ F21.4) 

DO 280 1=1/4 
AKT(P=<T(I) 

PC»100.*AKT( I ) /E ( 1/ 1 ) 

280 PRINT 230/ PC/ I 



230 FeRM*TV/F58*3# 11H PERCENT IN/I2*8H 00 BAND) 
RETURN 
END 



SUBROUTINE MATINIV ( A/ N, B , M, PETERM, NMAX ) 

DIMENSION At 4 , 4 ) , 3 ( 4 , 4 ) , P I V8T < 4 ) , INDEX! 1 *) 
DETERM=1 .0 
D8 20 1 = 1, N 
PIVOT! I ) =0,0 

20 INDEX! I 1=0.0 
D0 550 I = 1 , N 
AMAX=0.0 

DO 105 J=1,N 
IE iPIVOKj)) 105/21,105 

21 DO ID' -1 < = 1,N 

IF ( E I V 0 T ( K ! ) 100/22/100 

22 TEMP= A5S! A! J/K) ) 

IF (TEMP-AMAX) 100/23/23 

23 I R 0 '.•; = J 

I C8LUM=K 
AMAX=TEMP 
100 CONTINUE 
105 CONTINUE 

INDEX! I ) =4006* IR0W+ ICOLUM 
J= I ROW 

AMAX=A( J/ ICOLUM) 

DETI'RXIsAMAX*DETERM 
IF ( DETERM ) 24,600/24 

24 PIVOT! ICOLUM)=AMAX 

IF ( IRON- ICOLUM) 26/260,26 

26 DETERm*-DETERM 
DS 2C0 6=1, N 

SWA p =A(j,K) 

A ( J , < > A! ICOLUM, K) 

A! IC.5'_UM/K)aSWAP 
200 CONTINUE 

IF <M) 260,260,27 

27 DO 250 K = 1 , M 
SWAP* 3! J,K) 

B(J,K)=5(IC q LUM,K) 

B ( I C8LUM, K ) =SWAP 

250 CONTINUE 
260 K= ICOLUM 

A! I C8LUM, K ) a 1 • 0 
08 350 K = 1 , N 

A! ICOLUM, K)=A( IC8LUN,K)/AMax 
350 CONTINUE 

IF (M! 380,380,28 

28 DS 370 K*1,m 

B ( lCSLUM#K)aBt I C9LUM, K ) /AMAX 
370 CONTINUE 
380 D8 550 J = 1 , N 

IF (J-IC8LUM) 20,550,29 

29 T«A( J, ICSlUM) 

A! J, iCVLUM) =0.0 
DO 4 50 ,r = l,N 




At J,K)=A.< J,'<)-A< ICOLUM,K)»T 
450 C8NTINUE 

IF (M) 550,550.31 

31 00 500 K = 1,M 

B( J,K)=Q( IC0LUM,K)*T 

500 CONTINUE 
550 CONTINUE 
600 OB 710 1=1, N 
1 1 »N+ 1*1 

K= INDEXt I 1 J/A096 
IC8LUK= INDEXt II )-• 096*K 
IF (K-ICOLUN) 32,710,32 

32 DO 705 J = 1 , N 
SWAP- A( J,K) 

A(J,K)=\(J, I COLUM) 

At J, IC0LUM)=SWAP 

705 CONTINUE 
710 CONTINUE 
RETURN 
END 


o 





USERS MANUAL 


PROGRAM T e FIND BrST-FlJ EMPIRICAL DAMPING LAW THROUGH LEAST SQUARES METHOD 
READ (315) INDEX, 12, 10 

INDEX = 1 - FIND DAMPING LAW OF FORM D » C*X1»»A1 

- 2 D = C»X1»*A1*X2*#A2 

= 4 D = Al4«A2*Xl + A3*Xl**2 

= 5 (12 =7' D = C*m*»a1* »»A2*X**A3 

= 5 (12 =7) D = C<M**A1» ** A2* ( X/L ) ** A3 

12 .. 0 - xl s T 

= 1 - XI = X USED WITH INDEX = 1 

= 2 - XI = 

= 3 - Xi = M 

= 4 - XI = M, X2 * X 

= 5 - XI » M, X2 = USED WITH INDEX » 2 

= 6 - XI = X, X2 = 

= 7 - Xl = T, X2 = X/L (USED WITH INDEX = 2) OR X3 = X/L (INDEX » 5) 

IB = 0 - PRINT EACH DATA POINT 
= 1 - MINIMAL PRINT 

1 READ ( 13A6 » A 1, II) ( ARRAY ( I ) , 1 = 1,14), ISTOP 

ARRAY = TITLE OF DATA DECK 
Op ISTBP = 9 - STOP 

— . = 0 - CBNT’NUE 

2 READ < ID, I ID) 11,13 


= 

0 

- GO TO 1 

FREQ 

XN 

gen. mass 

F0RCE 

XE 

XO 

r 

1 

- INPUT UNITS 

• HZ 

H 

KC»-S**2/M 

KGF 

XE/XN 


s 

2 


HZ 

G/L5 

LB"S**2/ IN 

L3 

g/ls 


= 

3 


HZ 

M 

KG 

N 

XE/XN 


= 

4 


HZ 


N-M 

N-M 



= 

5 


HZ 

DEG 

K3-M**2 

N-M 

XE/XN 


= 

6 


HZ 

RAD 

LB-S**2/IN 

IN-KIP 

RAD/ IN-LB 


s 

7 


HZ 

g/lb 

L8“S**2/ IN 

LB 

g/lb 

IN 

13 

! = 

1 - 0 CALCULATED FR9M 0 

= pi*xe*f 





13 = 2 - D K 4»P I »ZETA*T 


READ ( E 10 . 3 ) CLGTH 

cloth = characteristic length of vehicle in meters 

3 READ (7E10.3) FREQ/ XN, GM, F, XE / X0» ZETA 


FREQ = FREQUENCY 

xn = d,;.flectibn at station used for normalization of gen, mass 



GM ^GENERALIZED MASS OR KINETIC ENERGY FOR II = A) 

F = EXCITATION FSRCF. (OR DISSIPATED ENERGY FOR II = A) 
XE = DEFLECTION AT EXCITATJ8N STATION 
XO = RADIUS AT NORMALIZATION STATION FOR 11=7 
ZETA = MODAL DAMPING COEFFICIENT (OR RADIUS (INCHES) AT 
FOR 11 = 7) 

IF (/FREQ/.LE-1.E-10) GO TO 2 

GO TO 3 • 


sample data deck 


1 0 

•313 

.75 

0 


TESTS NO. 11 

+ 13, SATURN UPPER 

STAGES 

bending tests 

3 

2 




33*858 





A. 33 

•1301 -2 

1 C-555 . 

3229. 

• 1733 

9.66 

.0335 -2 

160A8. 

3481 . 

.2341 

10.00 

.0290 -2 

19A21. 

3346. 

.3137 

10.58 

.0228 -2 

371 14. 

3449. 

.4005 

1 A . 80 

.0190 -2 

56828. 

3305. 

.6826 

19.01 

.0151 -3 

1053120. 

3412. 

5.169 

•23.1.2 

•Clog -3 

17980500. 

2763. 

16.13S 

A. DO 

• 1985 -2 

1023 A . 

3371 • 

.1764 

9.39 

• 361 t 3 

23731. 

3609* 

.2752 

1C. 90 

2 31 -3 

30539. 

3395. 

.2714 

13.33 

• 99 -4 

2284183. 

3325* 

.6684 

1 A . A 2 

• i 59 -3 

186813. 

35C2. 

• 7323 

16. <*0 • 

• 33 -4 

7 0 1 A A 3 . - 

3574. 

1 .6274 

22.97 

.19 -4 

8261987. 

4558* 

29.333 

5.20 

•2232 -2 

11824. 

3705* 

.6819 

9.63 

.393 -3 

1 3646 . 

3774. 

• 2203 

10.59 

.218 -3 

21531. 

34 38. 

• 3998 

10.91 

.163 -3 

22986. 

3736. 

.4452 

15.00 

•200 -3 

22630. 

3373. 

.5734 

2A.p6 

.570 -5 

82876450. 

4111* 

44.2231 

26.82 

.104 -4 

28463950. 

5444. 

71.8714 

5. 0 9 

.831 *3 

12523. 

3346. 

.9098 

9 . A A 

.546 - 3 

10630* 

3415* 

• 2459 

10.95 

• 3 5 7 - 3 

14796. 

3546. 

• 2593 

1 3. A9 

.151 -3 

554562. 

3732. 

.4576 

15.20 

• 158 -3 

21740. 

3907. 

• 3003 

3 

2 




33*853 





2. 0 9 

• i 0 8 7 -2 

973344. 

3q22 • 

• 4120 

7.81 

•9059 -2 

50200* 

5199. 

• 3559 

10*09 

.45 -5 

7571673. 

5466. 

4.0316 ’ 


2 


EXCITATION STATION 


*0089 

4 


1 

.0118 

4 


2 

.0123 

4 


3 

.0114 

4 


4 

.0130 

4 


5 


4 


6 


4 


7 

.0194 

4 


15 

.0109 

4 


16 

.0076 

4 


17 


4 


18 

*0095 

4 


19 


4 


20 


4 


21 

*0137 

4 


8 

.0094 

4 


9 


4 


10 

• 0115 

4 


11 

♦ 0108 

4 


12 


4 


13 


4 


14 

*0144 

4 


22 

•009 

4 


23 

*0114 

4 


24 

♦ 0100 

4 


25 

• 0109 

4 


26 


•o' 10 

3 - 

1 

• 0104 

3 - 

2 

• 0086 

3 - 

3 



* 


10.82 

' .141 -3 

158658* 

5482. 

.4760 


• 0163 

3 

. 

4 

14.31 

.140 -3 

28820* 

4478. 

.6195 


.OHO 

3 

• 

5 

1 7 * C2 

.14 -4 

1591757. 

4369. 

5-3168 



3 

- 

6 

2.08 

.864 .3 

1193429. 

3345. 

.3601 


.0180 

3 

- 

11 

7.63 

•3034 -2 

47160- 

3536. 

• 2568 


• 0084 

3 

- 

12 

10.38 

• 24 -4 

21984770. 

3782. 

5.0237 


•0072 

3 


13 

1 1 • D 1 

.126 - 3 

457090- 

3533. 

.3706 


• 0102 

3 

- 

14 

14.60 

.37 - 4 

859699. 

34 76* 

1.2591 



3 

- 

15 

1 7 . 22 

• 26 -4 

776740. 

3452* 

2*1896 



3 

- 

16 

8 . i 6 

•1645 -2 

21164. 

2327. 

• 2791 

2 

• .0098 

3 

- 

7 

1 H «q7 

.152 -3 

'9634. 

2330* 

.3545 


• 0085 

3 

- 

8 

S .06 

.1849 -2 

21705. 

2303. 

.2260 


• 0108 

3 

• 

17 

13.08 

« 4 4 - 4 

1136601 . 

2463. 

.7624 


• 0108 

3 

- 

18 

8.^9 

•1593 -2 

13517. 

1988. 

• 4762 

3 

• 0103 

3 

- 

9 

1**27 

• !.q7 -3 

15436. 

2277. 

.3594 


•OHO 

3 

- 

10 

3.31 

.2166 -2 

12720* 

201'* • 

.3786 


• 0099 

3 

- 

19 

13.60 

• 7o5 -4 

455873. 

2499. 

• 3541 



3 

- 

20 


TEST N0. 3 

SATURN v - 

DTV CONFIG. 

I T0RSI0NAL 

(1967) 



6 C 

2 







17-3 

5-029 








17-0 

5 . 660 

6*092 -5 

9*40 

+ 6 

219*12 

2.5 

-11 

.01 

17- 1 

6.362 

8*469 -5 

3 9 0 

+ 6 

227.04 

2.6 

-11 

•008 

17- 2 

8-117 

3*192 -5 

2*5 

+ 7 

20^-60 

2.4 

-11 

•016 

17- 3 

8.839 

1*436 -4 

3*88 

+ 6 

182.16 

4.5 

-11 

•01 

17- 4 

11.326 

6*703 -5 

2*7 

+ 6 

?4o»24 

8*0 

-11 

•02 

17- 5 

5.691 

6 • £7 1 -5 

9*1 

+ 6 

201 * 

3o • 

-12 

• 009 

17- 6 

■ 6.; 94 

6*432 -5 

6*8 

+ 6 

201 * 

30* 

-12 

• 009 

17- 7 

8.113 

3*293 -5 

3-1 

+ 7 

1 94 . 

40* 

-12 

• OH 

17- 8 

8.395 

1 .443 - 4 

4*0. 

+ 6 

205. 

1*C» 

-1? 

• Oil 

17- 9 

11.438 

A * 6 4 0 “5 

3.4 

+ 6 

238. 

50- 

-12 

♦ 018 

17-10 


9 



f ^endjob. - 

; ^assign g=.;t">si=cr,rq=kti#l?=lp, 

E* I ;.C v Ti- 
A"ORTtjA\ ;.G* US' • 


1 Dl-Er.3I 3N' '■•■>! t25] , 3MC253* TA3LEQI2S], T ABLEW T25] 

2 READ 25/ NM3DES/ IT^AX/ NQ, NX/ TBL 

2 23 rrifs^AT C 4I5/F10.3] 

4 ‘.'>9 27- ! ” J / NK«DES 

h 27 Read 35/ I]/ Gft C 13 

6 3: r * t r tors -33 


7 3 f. r O' 5 - AT tcF! 0*31 

' ~ E A D 33 / C/E.N- 

9 READ 30 / f. TABLED Cl] / I = 1 /NO] 

1C READ 30/ CTADLEWCI3/ I » 1/NW3 

i i print a o 

1? 40 F ftR- AT C1H1/P8V/4HFRE0/10X/9H3ENL MASS/ 7X/9HF0RC AMPL/ 7X 

13 * 5-FORC FRK3/7X/9HPESP aMPL/ 7X» IOhRESP PHASE//] 

14 DO 73 I = 1/ NN f *DES 

;5 DO 75 J = 1,1"/.' 

16 - TA2lE'/.’CJ3 

-.7 DO 75 :< s 1/ND 

i c j *v = T A 5 L E D C K 3 

l'“ C * 1 • 


I N r ' = 0 


r>,,.. T tt-> i .v, Reprodu 

■ PI .T -.l, . bes* av« 

1 '•! MeI J— — — — 

Cam, f r /K C I ] # v: , cmci]/ x, en, C/ on, t, phi] 

0 = X- T 

EPP * CQ -X]/X 

iF C A3S CERk] - TOL] 150/150/120 
' If [\ - 1 T f- A X ] 100/200/200 
. PRIM 160, >:N t J ] , SMI], ON/ W, 0/ PH» 

: F f;;>- AT C1PE36.4/5E16.4] 
r- pi t o 75 

. PRINT PIC, ITNAX/TOL/O/X 

. r 0' ,N 'A T C//16X/26HS0MJTISM DIDNT CONVERGE IN/I3/26H STEPS WITHIN TO 
»L£ D A\CE SF/F9.6, 7H XCN] = , 1PE10.3/ 9H XCN-1] »/E10.3] 

ST a P 

j CSNTINUF 



PROGRAM TO FIND RESFONSE OF NON LINEARLY 
DAMPED STRUCTURES 



= 

i 

SvWjukz f c<:\,w,sy.. 

= 

p 

'■ ?*•*•/. 

- 

■ 3 

T; = »-. ' . * a \‘ _ * W 2 

= 

»4 

Tl = T A « T A 

= 

c, 

T : 3 = 3" #30 

= 

£ 

TC = *?*TQ/2« 

= 

7 

TL = C7 1^/3.1^159 

= 

* S 

T F = T->TE 

= 

* 5*1 

T 3 = T r * * E *<’ 

= • 

10 

to = tg*tg 

= 

I 1 

tp = ro^ir 

= 

: 2 . 

= G Q P T C T 1 + T 2 3 

- 

* *j 

r = o - .'jk/ceno^ 



Apt's = K«T3/TA, 

= 

1 5 

= ATA.\‘CA:<3] 


16 

T 3 = 2. » CE *T2/3 

= 

: 7 

T 4 = DE.\'J?** r :E \‘' v *Dr;\ ; aM 


; 3 

r P = l . + CW3/T4 

= 

*. .T‘ 

T = r/rp 

= 

he 


= 

cl 

Eon 



P ?G3 A " ALL c A T 1 0 M 


o:o?6 

r 

03030 W? 

OUM'^Y 

w 

00032 

TA 

nUMI'Y 


00034 Tl 

00036 

TQ 

DU^Y 

GM 

0 J ^ ’ Y 

n 

CCQ*0 TC 

OCOAE 

TE 

DU^vy 

C 

coo.-* 

TP 

CCC*6 TG 

DU^y.y 

E\ 

00050 

TO 


T c 

00054 D£r%6M 

DL'^:‘'Y 

ON 

00056 

ARi 

V, v ‘ Y 

r- I 

OCCO'.O T3 

0006? 

T4 

OCOM 

FP 


°* POMVV T 

Sv8 p S?G'cA''S h L j J IP £ 0 

SCR7 ATAN 

ThE ElC 




USERS MANUAL 


PRQGRaM tS FIND RESPONSE 8F A MeNLINEARLY DaMPeD SYSTEM T9 A PERIODIC EXCITATION 
READ <4I5iF10.3> NMODES/ ITMAX, NO, NW, T8L 

NMODES = NO. OF MODAL EQUATIONS TO SOLVE 

itmax = maximum number of iterations for newtons method 

nq = no. of values of q umpl. of forcing fn.) for each forcing freq. 

NW = NO. OF VALU-S OF W (FREQ. OF FORCING FN.) 

tol - eR r o p tolerance requiring another newton iteration 

do 1 I = I/.;M3DES 

1 READ (2F8.3) ( W N ( I ) / GM(I)/ I = 1, NMODES) 

w.NU) = modal frequency of the i-th mode 

GM ( I > = Gen. MASS OF I-TH MODE 
READ < 2F8-3) C,EN ' 

C = COEFFICIENT OF EMPIRICAL DAMPING LAW 
EN = EXPONENT OF EMPIRICAL DAMPING LAW 

READ (8F10.3) ( TABLED! I ) / I = 1,NQ) 

OQ READ (8F10.3) (TASLEW(I)/ I = 1/NW) 

^ TAGLEQ(I) = VALUES OF Od) 

t a.elew ( i ) = values of wm (rad/seo 
sample data deck 

4 10 5 5 .001 

6.6266 8558*9 

11.0207 3121.6 

16*6826 ■ 4791 »6 

23.2190 10493.6 

♦313 »S 

1.-6 l *-4 1.-2 1 . 100 . 

6. 12* 18. 24. 30. 

ALL SATURN/ 747 DATA 

1 



aFNDjrS. 

* •' ss i g : - : n*y.T0/s:«cR>e9»MTi#L3«LP* 

A\r V. I :,D '-IT 1 • 

aFqRtRA'.' ['3,L':, 

* , i calculate &f equivalent linear Damping coefficient 

= 2 Z 

= 3 Dr-£\S!9N. TABLEC20] , TITLEC143 

- m P E A D 1 0 , C. , E 

-- ' 5 Ppr>iT 2,C,E 

= 62 F fj"N at C!M1/43X/25hDAMPING LAW CONSTANTS C =F7.3, 4H E =,F7.3] 

= 7 HE.-D 5, N A ’- 3 

= 3 5 FoT"AT C t C 1 5 3 

= 9 read ir... uasle tn, i = i,nampi 

= 10 10 r 0 r;v AT c ■ OF 3*31 

* 11 12 Pf'-D 5.NM3DES 

* 12 IF C- ••MODE S3 15,300# 15 

= 1 3 15 °EAD 20# CTITLEU3# I = 1#143 

= 14 2;; F 0 1 . AT C13-‘6»A23 

= 15 PRINT 20, 'TITLE C 1 3 # I = 1,143 

* 16 30 Fp'VAT ClriC>?9X, 13 a6# a2//2QX#4HM9DE» 5x,4hFREQ, 3X,9hGENL MASS 

= 1 7 * 7 >: » 9MA y PL ! TL'DE# 11X, 1HT, 15X, 1HD, 14X/4HZETA/ 26X# 9HIRAD/SEC3. 

= 15 * ?X, 4MCK33, 13X, 3H[m3, 12X, 5H[M-M3, 1 IX, 5H [N-M3 /3 

= 19 Dr; 75 i = 1 , N v 0DES 

= 20 R E O 50, W#GM 

= £1 50 F?p-iat Cl' El 0*3] 

= £ 2 T a • 5 * 3 ' A * a * W 

= £3 D“ 75 J = 1 # \AMP 

= 24 T i = T*TASLEU3*TA3LECJD 

= ?5 D = C » r 1 *-» E 

= 26 7. = D/12.56636/T1 

= 27 IF U-13 70,60,70 

00 = -,;5 60 r«?',T 65, I , W, DM, T A3LE C J3 , T1 , D, Z 

« 29 OeTfi 75 

* 30 '63 '' 6- : ';A r n23,r 12.3, 1P5E16-..4] 

* 21'. 7 ": PRINT 72, TA3LFCJ3 ,T1,0,Z 

= 22 7-2 FC-‘’"A7 C 1PF&7 ♦ 4 , GEl 6 • 43 

» 33 75 CO' T1 NOE 

« 34 GOTO 12 

= 35 30C ST^P 

* 36 End 


PROGRAM ALLOCATION 



oocoo 

T-'-BLE 

C0C50 

TITLE 

00104 

NAM P 

001C5 I 

v " ' 

00106 

N 'NODES 

00107 

J 

00110 

C 

00112 E 


00114 

V. 

CCU 6 

G M 

00120 

T 

00122 T 1 


CCl£4 

D 

C C 1 2 6 

1 






the end 
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PROGRAM TO CALCULATE EQUIVALENT 
LINEAR DAMPING COEFFICIENT \ 



USER'S MANUAL. 


program to calculate equivalent linear damping cceeficient from empirical 

DAMPING LAW D = C-»T**E 
READ ( 2E8 • 3 ) C,E 

C = C9EEF I C I ENT OF EMPIRICAL DAMPING LAW 
E = EXPONENT OF EMPIRICAL DAMPING LAW 

READ (15) N AMP 

name = no. of modal amplitudes per mode for which damping coeff. found 
read (IDES. 3) (TABlE<I)» I = 1/NAmP) 

TABLE! I) = VALUES OF MODAL AMPLITUDES 
1 READ (15) NMBDES 

NMODES * NO. OF MODES FOR WHICH DAMPING COEFF. TO BE FOUND 
READ (13A6/A2) ( T I TLE ( P > I ■ 1,14) 

title(I) = identification of structure examined 

DO 2 I = 1/ NMODES 
00 2 read (2E13.3) W.GM 

QQ W = FREQUENCY (RAD/SEC) of I-TH mode 

GM - GEN. MASS (KG) OF I-TH MODE 

GO TO 1- 

I F {NKRDES. £3.0 ) STOP 

SAMPLE DATA DECK 

.313 *8 

5 

1,-6 5 • *6 1.-5 5.-5 1,-4 

4 

SA 500D TIME POINT 2 ( DTV ) 

6.6266 8598*9 

11.0207 3121.6 

16-6826 4791.6 

23.2190 1C493.6 



